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Q_l' ABSTRACT 

A low mass-to-light ratio for the stellar component of spiral galaxies (M/L < 1 in 
the I-band) is advocated by various dynamical arguments and by recent cosmological 
I/*} ' simulations of the formation of these systems. We discuss this possibility by means of 

chemo-photometric models for galactic discs, adopting different Initial Mass Functions. 
We show that a number of "bottom-light" IMFs (namely, with less mass locked in 
low-mass stars than the standard Salpeter IMF), suggested independently in recent 
literature, do imply M/L ratios as low as mentioned above, at least for late type spirals 
(Sbc/Sc). This conclusion still holds when the bulge contribution to mass and light 
is included. We also predict the typical stellar M/L ratio, and correspondingly the 
zero-point of the Tully-Fisher relation, to vary considerably with Hubble type (about 
^5 ■ 0.5-0.7 mag in the red bands, from Sa to Sc type). 

| For some of the bottom-light IMFs considered, the efficiency of metal production 

. tends to exceed what is typically estimated for spiral galaxies. Suitable tuning of the 

IMF mass limits, post-supernova fallback of metals onto black holes or metal outflows 
must then be invoked, to reproduce the observed chemical properties of disc galaxies. 
In the appendix we provide M/L — colour relations to estimate the stellar M/L 
O ' ratio of a galaxy on the base of its colours, for several IMFs. 

U '• 

r/y , Key words: Galaxies: spirals; galaxies: chemical and photometric evolution; stars: 
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1 INTRODUCTION Navarro & Steinmetz 2001; Abadi et al. 2003; Salucci 2001; 

Salucci & Borriello 2001; Bottema 2002; de Blok & Bosma 

„ f 4.1. c 4- f j- 2002; de Blok et al. 2001a,b; Marchesini et al. 2002). 

Present cosmological simulations ol the formation of disc ' 

, . r, t^. , i-i ■ The suggestion of a low stellar M/L ratio is nevertheless 

galaxies withm the Cold Dark Matter hierarchical scenario, °° . ' 



reinforced by a simple estimate of the stellar mass of our 
Galaxy, based on the local stellar density and on a plausible 



suffer from the angular momentum problem (Navarro, Frenk 

& White 1995). Suggested ways out include early energy 
. . ,. c 4. n C j u i /c t r< i i density profile tor the Disc Sommer-Larsen & Dolgov 2001). 

injection from steffar feed-back (Sommer-Larsen, Gefato J ^ y ° ' 

& Vedel 1999; Sommer-Larsen, Gotz & Portinari 2003; The Stellar maSS ° btai T d in this way f ° r the Milky Way is 

Thacker & Couchman 2001) or Warm Dark Matter cosmol- of the ° rder ° f 5 X 10 ' and if ° ur Galaxy is to lie 

/0 t c t~. i nnm no i_ M u j , unn on the observed TF relation together with other spirals, its 

ogy (Sommer-Larsen & Dolgov 2001). Such N-body+SPH ° p '. 

, .. , . .„ , M/Lj ratio must be around 0.75 (see lig. 1). Considering 

simulations not only improve significantly on the angular ' v ° \ 

, n , . , » I, . ,f all uncertainties, Sommer-Larsen & Dolgov (2001) estimate 

momentum problem, but also compare successtulfy to the 3 ; ' 

u j m. ii T-i- i. jtm i4- j j 41 for the total baryonic matter in the Milky Way (inclusive of 

observed lully-1 isfier ( 1 r ) refation, provided the mass-to- ■' j j \ 

v i 4 /tvi it \ 4- r 41 4 ii 4-i ic< stars, gas and stellar remnants) M/L/ = 0.87 ± 0.57. 

light (M/L) ratio of the steffar component is fow Isommer- ' ° 11 

t o t~. i nnm o t 41 nnnoN • 41 t Similar conclusions can be drawn from the Sb galaxy 

Larsen & Dolgov 2001; Sommer-Larsen et al. 2003): m the 1- , T „„ r . 

u j4i.-r jai/t no t \ it ■ NGC 2841: estimates of its disc mass range between 

band, the interred M/L/ ~ 0.8 (rig. 1). However, m cosmo- , 1n , ,„ . ° 

,/ a t \ \ ■ A/r/T 4- x 10 10 M Q and 11.7 x 10 10 M * and its circu- 

logical simulations the need of a low baryonic M/L ratio may M u u 

be partly due to the too dense cores of dark matter haloes: 

in simulated galaxies, dark matter is more concentrated and 

dynamically dominant in the central regions, than estimated * We indicate with h = h 100 , with h 50 and with h 75 the Hubble 

from rotation curves (Navarro & Steinmetz 2000a, b; Eke, constant in units of 100, 50 and 75 km sec -1 Mpc -1 
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lar velocity is 317 km sec - (Kent 1987; Benson et al. 2000 
and references therein), which also points to a M/L/ < 0.7 
(Fig. 1). Regarding the zero-point of the Tully-Fisher rela- 
tion, however, some uncertainties exist and are discussed in 
Appendix A. 

An independent indication of a low M/L ratio for the 
baryonic component of spiral discs relies on stability argu- 
ments. Bar instability of isolated discs was one of the reasons 
to introduce massive dark halos historically (Ostriker & Pee- 
bles 1973); the argument still holds, in the sense that if too 
much mass is stored in the luminous disc component, strong 
bars should form. Based on this, Efstathiou, Lake & Negro- 
ponte (1982) put an upper limit of M/Ls < 1.5ft to the disc 
M/L ratio, that is M/L B < 1 for h = 0.65. 

The M/L ratio of the stellar component is also relevant 
to the issue as to whether discs are maximal or sub-maximal, 
or to what extent the luminous, baryonic mass dominates 
the dynamics and the rotation curve in the inner regions of 
the disc. Even under the hypothesis of maximal discs, Bell & 
de Jong (2001) point out that lower M/L ratios for the stellar 
component are implied, than those predicted with a Salpeter 
IMF. And discs may well be sub-maximal (Bottema 1993, 
1997; Courteau & Rix 1999), although this is still much de- 
bated (Ratnam & Salucci 2000; Fuchs 2002, 2003a,b; Masset 
& Bureau 2003; Kranz, Slyz & Rix 2003). For his favoured 
sub-maximal disc case, with the luminous component con- 
tributing 63% of the rotation curve at maximum, Bottema 
(2002) suggests a stellar M/L/ ~ 0.82 for spirals in the Ursa 
Major cluster. For the Milky Way itself, a variety of mass 
models is compatible with available kinematic data (Binney 
& Merrifield 1998; Dehnen & Binney 1998), so that on this 
base no clear discrimination is possible between the maximal 
or sub-maximal case; though recent combined results on lu- 
minosity and mass distribution, gas flows and microlensing 
seem to favour a maximal disc (Gerhard 2002). 

Two recent, detailed studies of individual spirals also 
point toward low M/L ratios. A dynamical analysis of the 
isolated Sc galaxy NGC 4414 indicates an intrinsic stellar 
M/L<~1 in the B, V and I band (Vallejo, Braine & Baudry 
2002). For the Sab spiral 2237+0305, Huchra's lens, Trott & 
Webster (2002) combine dynamical and gravitational lensing 
constraints to derive M/L/ ~1.1 for the disc. 

In the light of these results, the obvious question arises 
whether M/L ratios so low are plausible from the point 
of view of the stellar content and of the photometric evo- 
lution of galaxies; this relies heavily on the stellar Initial 
Mass Function (IMF) of stars in disc galaxies. A standard 
Salpeter IMF, extended over the typical stellar mass range 
[0.1-100] Mq certainly yields much higher M/L ratios than 
those mentioned above (see §5). On the other hand, there 
is by now a general consensus that a single-slope power 
law IMF over the whole stellar mass range is not adequate 
for the local nor for other environments. The slope of the 
IMF appears in fact to become shallower below 1 Mq, in 
the local environment (Miller & Scalo 1979; Scalo 1986; 
Kroupa, Tout & Gilmorc 1993; Reid & Gizis 1997; Gould 
et al. 1997, 1998; Chabrier 2001) as well as in globular clus- 
ters and in the Galactic Bulge (Piotto & Zoccali 1999; Zoc- 
cali et al. 2000). Besides, a turn-over in the regime of brown 
dwarfs is favoured by recent results (Chabrier 2002). The 
observed IMF, that is, appears to be "bottom-light" with 
respect to a single power-law, Salpeter-like IMF; for recent 
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Figure 1. Straight lines: location of the observed Tully— Fisher 
relation for Sbc-Sc disc galaxies by Giovanelli et al. (1997a,b; 
h=0.65 adopted here) in the stellar mass vs. circular velocity 
plane, for different assumptions about the stellar M/Lj ratio. We 
show also the location of simulated disc galaxies by Sommer- 
Larsen et al. (2003, triangles), and of two observed galaxies: the 
Milky Way and NGC 2841 (the latter rescaled to h = 0.65). 

reviews, see Kroupa (2002); Chabrier (2003). From the theo- 
retical point of view, the flattening and turn-over of the IMF 
at low masses may be related to the Jeans scale (thermal 
support) and to the scale of magnetic support against grav- 
itational collapse (Larson 1998; Padoan & Nordlund 2002). 

In this paper we investigate what kind of IMF is re- 
quired to obtain a low M/L ratio (say, ~0.8 in the I-band) in 
a typical Sbc-Sc galaxy. We consider different IMFs among 
those suggested in literature from a variety of independent 
studies. For each of these IMFs, we calculate chemical evo- 
lution models with infall, metallicity gradients, and star for- 
mation histories representative of late-type spiral discs; then 
we compute the corresponding broad-band luminosities by 
convolving suitable Single Stellar Population (SSP) models 
with the star formation and metal enrichment histories. 

The paper is organized as follows. In Section 2 we dis- 
cuss which Hubble types and photometric passbands are 
most relevant for our present analysis of the M/L ratio in 
galactic discs. In Section 3 we present the various IMFs con- 
sidered in our chemo-photometric models and in Section 4 
the corresponding sets of SSPs. In Section 5 we present sim- 
ple exponential models with fixed metallicity, to give a first 
overview of the problem and to discuss the solidity of the 
predicted M/L ratios with respect to the adopted stellar 
models and spectral library. In Section 6 we introduce our 
full chemo-photometric models for disc galaxies, and the 
corresponding calibration on observational constraints. In 
Section 7 we present our results on M/L ratios obtained from 
chemo-photometric models with different IMFs. In Section 8 
we add the contribution of the bulge to the global M/L ra- 
tio of galaxies. Finally, in Section 9 we outline summary 
and conclusions. In Appendix A we discuss observational 
uncertainties in the zero-point of the TF relation. In Ap- 
pendix B we provide M/L — colour relations to estimate the 
stellar M/L ratio of a galaxy on the base of its colours. 



© 2003 RAS, MNRAS 000, 000-000 



On the M/L ratio and the IMF in disc galaxies 3 



2 FOCUSING ON LATE HUBBLE SPIRAL 
TYPES 

The M/L ratio of the stellar component of a galaxy depends 
on two basic quantities: the IMF and the star formation 
history (SFH); the latter sets the relative weight of young, 
bright and blue stellar populations with respect to old, dim 
and red ones. The SFH depends on Hubble type, as clearly 
indicated by the variation of typical colours along the Hub- 
ble sequence (Roberts & Haynes 1994). Hence, the M/L ratio 
of the stellar mixture is also expected to depend on Hubble 
type. This effect is indeed observed in the empirical TF re- 
lation, where a systematic, negative offset in luminosity is 
observed when going from Sbc-Sc spirals to Sb to earlier 
types (Rubin et al. 1985; Giovanelli et al. 1997a); and the 
offset correlates with the B-R colour as expected from dif- 
ferent SFHs in different Hubble types (Kannappan, Fabri- 
cant & Franx 2002). Giovanelli et al. (1997a,b) correct for 
this morphological dependence, and provide a TF relation 
(the one shown in Fig. 1) typical for late spiral types. This 
allows us to focus on the M/L ratio of Sbc-Sc spirals. 

Having selected our fiducial Hubble type, we need to 
characterize it for the sake of the chemo-photometric mod- 
els to be developed. The historical, B-band based picture of 
the Hubble sequence as a sequence of bulge-to-disc (B/T>) 
ratios is not supported by recent surveys in redder bands; 
in the I or K bands, which are better mass tracers, B/T> 
correlates fairly weakly with spiral type, and with a huge 
scatter, thus proving to be a poor morphological indicator 
(de Jong 1996a; Sommer-Larsen et al. 2003 and references 
therein; see also §8). Kennicutt, Tamblyn & Congdon (1994, 
hereinafter KTC94) demonstrated that the sequence of spi- 
ral types is rather a sequence of different SFHs in the discs 
themselves, as traced by the so-called birthrate parameter 



SFR 

< SFR > 



(1) 



defined as the ratio between the present and past av- 
erage star formation rate (SFR). KTC94 estimated the 
fa-parameter both on the base of the equivalent width of 
the H a line, and from the ratio between the H a luminos- 
ity (tracer of the present SFR) and the V-band luminosity 
(tracer of the global underlying stellar mass in the disc) . The 
two methods are not completely independent, for EW(H Q ) is 
also a flux ratio between the H a line and the nearby red con- 
tinuum (corresponding to the R-band, roughly), but they 
give consistent results. From KTC94, the 6-parameter in- 
creases from b < 0.2 for Sa-Sab discs, to b ~0.3-0.4 for Sb 
discs, to b ~0.8-l for Sbc-Sc discs, while the trend breaks 
down toward later Hubble types, which display more irreg- 
ular and episodic star formation histories; see also Sommer- 
Larsen et al. (2003). 

Henceforth, in gross terms our chemo-photometric 
models will be taken to represent typical Sbc-Sc spirals 
when their SFHs correspond to a ^-parameter ~0.8-l. 



2.1 Focusing on the I— band 

In this paper we will mainly discuss the M/L ratio in the 
I-band, firstly because the observational evidence quoted in 
§1, including the fiducial TF, refers to the I-band. Secondly, 
red bands like the I-band are good mass tracers because they 




Figure 2. M/L ratio as a function of age for SSPs of differ- 
ent metallicity. Within this wide metallicity range (from ~ 10~ 2 
to 3.5 times solar), the M/L ratio in the I-band varies only 
by a factor of 2 (lower panel) as opposed to a factor of 5 
in the B-band (upper panel). 

These M/L ratios are computed from the mass of living stars on 
the SSP, not including the remnants of more massive, dead stars. 
The mass of remnants is however included in the exponential and 
chemo-photometric models in the remainder of the paper. 



are less sensitive than bluer bands to recent SF episodes, as 
well as to metallicity effects (Fig. 2). Thirdly, compared to 
redder, near infrared bands, the integrated I-band luminos- 
ity at intermediate and old ages is not so sensitive to the 
specific treatment of the Asymptotic Giant Branch (AGB) 
phase (see §5.2). Finally, red bands like the I-band are less 
affected by dust extinction and corresponding corrections. 
The empirical TF relation by Giovanelli et al., as well as 
other observational evidence quoted in the introduction, are 
already corrected for extinction. 



3 THE INITIAL MASS FUNCTION(S) 

In this section we discuss the IMFs we have considered 
in turn in our chemo-photometric models, together with 
their respective theoretical/observational background. We 
will present sets of models with three "power-law IMFs" 
and three "exponential IMFs", displayed in Fig. 3. In the 
following, M will indicate the stellar mass in solar units. 

(1) A Salpeter (1955) IMF, extended over the typical stel- 
lar mass-range [0.1-100] Mq: 

-1.35 



$(M) = C s M~ 



C s = 0.1716 



where C s is a normalization coefficient fixed so that the IMF 
is normalized to unit mass when integrated between the low 
and high stellar mass ends. 

Although the original Salpeter IMF was derived only for 
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Figure 3. Comparison between the 6 different IMFs consid- 
ered in this paper, all normalized to and integrated mass of 
1 Mq. Top panel: the three "power-law" IMFs (1—3), defined 
in the range [0.1-100] Mq. Bottom panel: the three "exponen- 
tial" IMFs (4-6), defined in the range [0.01-100] Mq; Larson (1) 
and Larson (2) refer to the original and "modified" Larson IMF, 
respectively; the Salpeter IMF is also repeated for comparison. 



stars between 0.4 and 10 Mq and, as mentioned in the in- 
troduction (§1), there is substantial evidence that a single 
power-law over the entire mass range is not adequate for the 
observed IMF, this functional form is still widely adopted in 
literature and it is useful to include it in the present study, 
for the sake of comparison to other IMFs. 

(2) The Kroupa (1998) IMF, derived for field stars in the 
Solar Neighbourhood and often adopted in chemical evolu- 
tion models for disc galaxies (e.g. Boissier & Prantzos 1999): 



$(M) 



Ckri = 0.480 



C krl M" - 5 
C kr M- 1/2 

CkrM- 1 - 7 



0.1 < M < 0.5 
0.5 < M < 1 
1 < M < 100 

Ckr = 0.295 



where Ckri and Ckr are fixed so as to guarantee normaliza- 
tion and continuity of the IMF over the range [0.1-100] Mq. 
This IMF is steeper than the Salpeter one in the high- 
mass range (M > 1 Mq) but it flattens out progressively at 
low masses (M < 1 Mq). The steep slope in the high-mass 
range was actually taken after Scalo (1986). In a more re- 
cent determination Kroupa (2001) finds instead a shallower 
slope of 1.3, close to the Salpeter value, but the steeper Scalo 
slope is recovered if unresolved binary systems are taken into 
account (see also Kroupa 2002). 

(3) The Kennicutt (1983) IMF, often advocated as ade- 
quate to reproduce the global properties of spiral galaxies 



(e.g. KTC94, Sommer-Larsen 1996), again with the typical 
mass range [0.1-100] Mq. 



$(M) 



Cfc83 = 0.328 



C kS 3 M 
C fe8 3 M 



-0.4 
-1.5 



O.K M < 1 
1 < M < 100 



This IMF is somewhat shallower in the high-mass range 
than the local one, as seems to be required to account for the 
observed H„ luminosities and equivalent widths of external 
disc galaxies. The flattening below 1 Mq is introduced after 
the observational evidence for the local IMF (see KTC94 for 
details). 

(4) The Larson (1998) IMF, over the mass range 
[0.01-100] Mq, that is down to the sub-stellar regime (which 
in this case has, however, a negligible contribution to the 
mass budget): 

*(M) = C L1 M- 1 - 35 exp(- Ml1 * 



M 



Cli = 0.317 



Mr 



0.3375 



This IMF recovers a Salpeter power-law at high 
masses, while at low masses the exponential cut-off 
produces a characteristic peak mass, or mass scale, 
M p = Mii/1.35 = 0.25 Mq in agreement with some obser- 
vational evidence for the Solar Neighbourhood. An exponen- 
tial cut-off in the low-mass range is favoured also by recent 
determinations of the local IMF down to the brown dwarf 
regime (Chabrier 2001, 2002; see below). 
This functional form of the IMF was originally introduced 
by Larson (1986) with the suggestion that the characteristic 
mass scale Mli may be larger for higher temperatures of 
the parent gas, as expected for the Jeans mass scale. This 
would imply a top-heavy IMF in the early galactic phases, 
at high redshifts and/or in low-metallicity environments. A 
suitable variation of the typical mass scale Mli could re- 
produce the systematic variations of M/L ratios in elliptical 
galaxies (Chiosi et al. 1998) and explain the chemical en- 
richment of the intra-cluster medium (Chiosi 2000; Moretti, 
Portinari & Chiosi 2003; Finoguenov, Burkert & Bohringer 
2003), or be responsible for the possible presence of dark 

massive objects in the Galactic halo (Chabrier 1999)T. 
In the present study however, we do not consider variations 
in the IMF with time or metallicity, and within each chemo- 
photometric model we assume the IMF to be constant. 
(5) A modified Larson IMF: 



$(M) = CL2M- 1 - 7 exp(- 
C L 2 = 0.4337 



1 



ML2 

M J 
M L2 = 0.425 



This IMF maintains the same functional form as the original 
Larson (1998) IMF, but with a steeper slope at high masses 
in better agreement with local determinations (Scalo 1986; 
Kroupa 2001). The mass scale Ml2 is chosen so as to main- 
tain the same peak mass M p — Mr,2/1.7 = 0.25 Mq as in 
the previous case. 

t Notice that the Macho— inspired halo IMF suggested by 
Chabrier (1999) and mentioned here, peaking at 1.7 Mq, dif- 
fers from the IMF derived by Chabrier (2001, 2002) for the Solar 
Vicinity, peaking around 0.1 Mq and used below in this paper. 



© 2003 RAS, MNRAS 000, 000-000 



On the M/L ratio and the IMF in disc galaxies 5 



(6) The IMF suggested by Chabrier (2001, his case IMF3): 



$(M) = Cc M~ 2:i exp 



_ (McV 
\ M J 



C c = 40.33 



M c = 716.4 



over the mass range [0.01-100] Mq. Chabrier (2001) adopts 
this functional form, a power-law with an exponential cut- 
off, after Larson (1986) and tunes its parameters to fit the 
local field IMF in the low-mass range (M < 1 Mq), showing 
also that this form of the IMF is valid down to the brown 
dwarf regime (Chabrier 2002) . In the range of massive stars 
the steep Scalo/Kroupa slope is recovered. 





Girardi/Castelli 


Lcjcunc 


bol 


4.720 


4.720 


U 


5.666 


5.576 


B 


5.448 


5.446 


V 


4.788 


4.807 
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4.413 


4.433 
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4.065 


4.080 


J 


3.640 


3.653 


H 


3.295 


3.307 


K 


3.263 


3.252 



Table 1. Solar magnitudes in the main Johnson-Cousins-Glass 
bands as predicted by the stellar libraries in use 



4 THE SINGLE STELLAR POPULATIONS 

A Single Stellar Population (SSP) is a coeval, chemically 
homogeneous assembly of stars with masses distributed ac- 
cording to the IMF; the integrated spectrum and luminosity 
Lssp(t, Z) of a SSP of age r and metallicity Z can be com- 
puted on the base of stellar isochrones, combined with a 
spectral library. SSPs are the "building blocks" of spectro- 
photometric models of more complex systems, like galaxies: 
by convolving SSPs with the SFH and the metal en- 
richment history Z(t) provided by chemical models, the cor- 
responding global luminosity in any given pass-band can be 
calculated as: 

L= [ 9(t)L gS p[T-t,Z(T-t)]dt (2) 
Jo 

where T is the present age of the galaxy. 

For the purpose of this paper, we computed SSPs 
for the IMFs discussed in §3, and for metallicities 
Z=10- 4 , 4 x 10~ 4 , 0.001, 0.004, 0.008, O.O19(=Z ), 0.04, 
0.07, 0.1. Our SSPs are based on the latest Padua 
isochrones (Girardi et al. 2000, 2002; Salasnich et al. 2000; 
http://pleiadi.pd.astro.it), recently extended to Z=0.1 by 
L. Girardi (private communication). These isochrones are 
constructed from stellar tracks with updated opacities, equa- 
tion of state and neutrino cooling, and with a moderate 
amount of convective overshoot; see the above mentioned pa- 
pers for all details. In this work we consider only isochrones 
with solar-scaled metallicity; a-enhanced compositions may 
be found for stars at low metallicities and/or in the inner 
disc regions, where the SFH is more rapid (see §6). Recently, 
stellar tracks and isochrones with a-enhanced composition 
have become available for Z >0.008 (Salasnich et al. 2000); 
for lower metallicities, suitable conversion formulae from the 
a-enhanced to the solar-scaled case can be applied (Salaris, 
ChiefR & Straniero 1993). We will consider the role of a- 
enhanced compositions in future models. 

For the sake of the present study, most notable in these 
isochrones is that the AGB phase includes up-to-date pre- 
scriptions for the Thermally Pulsing (TP) regime (Girardi 
& Bertelli 1998); the contribution of AGB stars is very im- 
portant for the evolution of the NIR flux and of visual- 
NIR colours (Girardi & Bertelli 1998; Schulz et al. 2002; 
Mouhcine & Langon 2002), and is potentially interesting also 
for the red bands considered in the present paper. The ef- 
fect of including full semi-analytical models of the TP- AGB 
phase in the isochrones and SSPs will be discussed in §5.2, 



by considering the isochrones by Marigo & Girardi (2001) 
as an alternative set. We stress that a consistent inclusion 
of all post-main sequence phases of stellar evolution in the 
computation of isochrones and SSPs is crucial for a cor- 
rect estimate of the integrated luminosity and of M/L ratios 
(Buzzoni 1999, 2002). 

The transformation from theoretical quantities (Ltoi, 
T e ff, g) into magnitudes and colours is performed through 
the spectral library compiled by Girardi et al. (2002), a suit- 
able extension of Kurucz' ATLAS9 library (Castelli, Gratton 
& Kurucz 1997). From the same library we derive for the so- 
lar model of the Padua tracks (log T e //=3.762, log g=4.432, 
L. Girardi, priv. comm.) the solar absolute magnitudes listed 
in Table 1, first column. To express integrated luminosities 
of photometric models in solar units, we consider in fact 
more consistent to refer to the theoretical solar magnitudes 
derived from the same spectral library, rather than to the 
observed ones or to the magnitudes given by more detailed 
solar spectrum models; differences are anyways well below 
0.1 mag in the bands of interest. Uncertainties related to 
the use of different libraries are briefly discussed in §5.3, by 
means of SSPs calculated with the Lejeune library. 

For further details on the calculation of SSPs, the reader 
is referred to Bressan, Chiosi & Fagotto (1994) and Tan- 
talo et al. (1996). For consistency with previous work, in 
the SSPs the adopted solar bolometric magnitude and the 
definition of the bolometric corrections are slightly different 
than those in Girardi et al. (2002). However, when allowance 
is made for the different definition of the zero points, the in- 
tegrated luminosities are identical in the two cases and the 
predicted M/L ratios are unaffected, as expected. 



5 A DISCUSSION OF SIMPLE MODELS 

For a preliminary study of the effect of different IMFs on the 

M/L ratio, we first resort to very simple models, adopting 

an exponentially decaying star formation rate 

_ t 
^(t) tx e t sf 

where different timescales tsf correspond to different val- 
ues of the 6-parameter (Eq. 1). The age of the model discs 
is assumed to be T=10 Gyr, as estimated for the bulk of 
the stellar population in the disc of the Milky Way (Carraro 
2000). Exponentially decaying SFRs with varying timescales 
have often been used as effective prescriptions to reproduce 
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Figure 4. M/L ratio in the I— band for simple exponential models with varying 6-parameter; different panels correspond to models with 
different IMFs, as indicated. The shaded area marks the range M/Lj=0.7— 1 favoured by observations of Sbc— Sc discs (b = 0.8 — 1). 



the Hubble sequence in spectre— photometric models (e.g. 
Larson & Tinsley 1978; Bruzual 1983; Guiderdoni & Rocca- 
Volmerange 1987; Bruzual & Chariot 1993; KTC94; Moller, 
Fritze-v.Alvensleben & Fricke 1997; Bell & de Jong 2000). 
Our preliminary "exponential" models do not include chem- 
ical evolution, but we can learn about the role of metallicity 
by adopting SSPs of different metallicities in Eq. 2. These 
simple models are also useful to assess uncertainties related 
to the detailed prescriptions for the TP-AGB phase and to 
the adopted stellar library (§5.2 and §5.3). 

Fig. 4 shows the stellar (i.e. stars+remnants) M/L ratio 
in the I-band for our exponential models as a function of 
the corresponding 6-parameter; different panels correspond 
to models with different IMFs. Metallicities between ~half 
solar and solar are considered representative of the bulk of 
stars in galactic discs. The shaded area indicates the range 
of observed M/L/=0.7-l we want to reach with our Sbc-Sc 
disc models (b = 0.8 — 1, see §2). 



It is immediately apparent how a Salpeter IMF yields 
a M/L ratio much higher than desired. The local Kroupa 
IMF has still a slightly high M/Li > 1 while all the other 
"bottom-light" IMFs generally reach M/L/ <1 for the mod- 
els corresponding to Sbc/Sc spirals. This overview promotes 
the development of more refined models for a detailed pre- 
diction of M/L ratios with bottom-light IMFs (§6). 

Fig. 5 shows for two examples, the Salpeter and the 
Kennicutt IMFs, the predicted M/L ratios of the global 
stellar population (i.e. stars and remnants, thick lines) com- 
pared to the M/L ratio of the sole population of living stars 
(thin lines). Neglecting the mass of stellar remnants (white 
dwarfs, neutron stars and black holes) induces an underesti- 
mate of the M/L ratios typically of 10%, which can increase 
up to 30% depending on the IMF and SFH (cf. models with 
low ^-parameter and Kennicutt IMF). This illustrate the 
importance of calculating self-consistently the contribution 
of stellar remnants when discussing the M/L ratio of the 
stellar/baryonic component of galaxies. 
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Figure 5. I-band M/L ratio for the global stellar population 
(stars+remnants, thick lines) compared to that of the sole living 
stars (thin lines). The mass of remnants typically contributes for 
a 10%, but it can reach a 30% for bottom-light IMFs, such as 
the Kennicutt IMF, and SFHs skewed to old populations (low 
b- parameters) . 



5.1 Dependence of the M/L ratio on Hubble type 

It is interesting to consider the variation of stellar M/L ratios 
with Hubble type, here parameterized by the 6-parameter 
indicative of the different star formation histories. Fig. 6 
shows the M/L ratio as a function of 6, normalized to the 
value corresponding to 6=1. The scale on the right axis in- 
dicates the corresponding shift in magnitude. The plot is for 
different IMFs and for solar metallicity; results for relative 
M/L ratios do not change considering the case Z = 0.008 
instead. In the I-band (top panel in Fig. 6) with respect to 
Sbc-Sc spirals we predict the M/L ratio to increase by a 
factor of 1.3—1.5 for Sb spirals (6 ~0.35) and of 1.7-2.2 for 
Sa-Sab spirals (6 ~0.1). This corresponds to a systematic 
magnitude offset of 0.3-0.4 mag and of 0.6-0.8 mag for Sb 
and for Sa-Sab spirals, respectively, as a result of the dif- 
ference in characteristic SFH. These offsets are somewhat 
reduced when bulges are included in the computation of the 
total M/L ratios of galaxies (§8.1), but they remain signif- 
icant. The offsets we predict are larger than the empirical 
ones found by Giovanelli et al. (1997a): 0.1 mag for Sb spi- 
rals and 0.32 mag for earlier types. However, the extent of 
systematic offsets in the TF relation is still a matter of de- 
bate: for instance, larger offsets in the nearby R-band are 
found by Kannappan et al. (2002). Their offset of 0.76 mag 
for Sa galaxies is in good agreement with our predictions in 
the R-band at 6 ~ 0.1 (Fig. 6, lower panel). 



Figure 6. Relative M/L ratio normalized to the value corre- 
sponding to b = 1 models. Top panel: I-band; bottom panel: ri- 
band. Results are plotted for different IMFs as indicated, for solar 
metallicity. 



5.2 The treatment of the AGB phase 

For SSPs of intermediate age, the luminosity (in particular 
the bolometric and NIR) is dominated by the contribution 
of stars in the AGB phase. Hence careful modelling of these 
phases of stellar evolution is very important when determin- 
ing M/L ratios, especially in the NIR. 

As underlined in §4, our SSPs are based on isochrones 
(Girardi et al. 2002) which explicitly include the TP- AGB 
phase with a synthetic algorithm based on the core growth, 
the core mass-luminosity relation and the initial-final core 
mass relation (Girardi & Bertelli 1998). Recently, Marigo & 
Girardi (2001) discussed self-consistency between chemical 
and photometric evolution in galactic models. They refined 
the calculation of the "stellar fuel" (amount of hydrogen 
or helium consumed by nuclear burning) adding the contri- 
bution of the chemical yields dredged up and lost by the 
stars, to the classic fuel estimate based only on the initial- 
final core mass relation. According to the Fuel Consump- 
tion Theorem (Renzini & Buzzoni 1986), this additional 
fuel corresponds to additional emitted luminosity (mostly 
in the NIR for AGB stars). Marigo & Girardi also provide 
isochrones including detailed semi-analytical models for the 
TP-AGB phase, improving upon the simple synthetic algo- 
rithm adopted by Girardi et al. 

In this Section, we check the effect of a different mod- 
elling of the TP-AGB phase on the predicted M/L ratios 
in galactic discs. We calculated a set of SSPs based on the 
isochrones by Marigo & Girardi (M-SSPs) to be compared 
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Figure 7. Comparison between the G-SSPs and the M-SSPs 
(thick and thin lines, respectively) in bolometric, I and K magni- 
tude. Detailed modelling of the TP-AGB phase (M-SSPs) leads 
to a larger overall luminosity contributed by AGB stars, mostly 
in the K and NIR bands. 

For the M-SSPs with Z = 0.008, notice the occurrence of a "red 
spike" around 1.5 Gyr, a feature discussed by Girardi & Bertelli 
(1998). This feature is very short lived and does not influence the 
bulk of emitted light integrated over time. 



to our reference SSPs based on the isochrones by Girardi 
et al. (G-SSPs). Fig. 7 compares the luminosity evolution of 
the G-SSPs (thick lines) to that of the M-SSPs (thin lines) 
for the three metallicities available for the M-SSPs; we adopt 
a total SSP mass of 1 Mq and the Salpeter IMF for the sake 
of illustration. With the M-SSPs, although the onset of the 
AGB contribution is more gradual (with a smoother phase 
transition due to the overluminosity effect, see Girardi & 
Bertelli 1998) , the overall luminosity produced in the course 
of evolution is larger, as expected from the larger "fuel" 
actually consumed. The effect is generally negligible in the 
I-band (and in all other optical bands) with the exception of 
very low metallicities, but it is highly significant in the NIR; 
in the K-band, the M-SSPs can be up to 0.6 mag brighter. 

Fig. 8 compares the M/L ratios of our simple exponen- 
tial models for discs, as predicted with the G-SSPs and with 
the M-SSPs. In the I-band the difference is minor, and neg- 
ligible for the relevant metallicities Z > 0.008, while in the 
K-band it becomes significant (up to 20-30%). The Salpeter 
IMF is adopted for the sake of illustration, but the situation 
is similar for all the other IMFs. 

In summary, the effect of detailed TP-AGB modelling 




Figure 8. M/L ratios for our exponential disc models, as ob- 
tained with the G-SSPs and the M-SSPs (thick and thin lines, 
respectively). 



is negligible in the I-band (and in optical bands in general) , 
making from this point of view the I-band luminosity supe- 
rior to NIR luminosities for mass estimates. 

Stellar and photometric evolution in the TP-AGB 
phase can be further altered by dusty circumstellar en- 
velopes (Bressan, Granato & Silva 1998; Mouhcine 2002; Pi- 
ovan, Tantalo & Chiosi 2003) or different opacities and mass 
loss rates in the carbon star phase (Marigo 2002; Marigo, Gi- 
rardi & Chiosi 2003). Again, the effects are seen mostly in 
the NIR colours but they should be negligible for the in- 
tegrated I-band luminosity, where the contribution of AGB 
stars is not so dominant. We conclude that estimates of M/L 
ratios in the I-band are robust with respect to uncertainties 
in the modelling of the TP-AGB phase. 



5.3 An alternative spectral library 

An alternative spectral library has been recently provided by 
Lejeune, Cuisinier & Buser (1997, 1998) and has been widely 
used in population synthesis studies (Boissier & Prantzos 
1999; Liu, Chariot & Graham 2000; Schulz et al. 2002). This 
is also an extension of the synthetic Kurucz' ATLAS9 li- 
brary, but with additional empirical corrections to reproduce 
the observed color-temperature relations of individual stars. 
However, Girardi et al. (2002) argue against a posteriori cal- 
ibrations of theoretical libraries, while Westera et al. (2002) 
find that it is not possible to attain a semi-empirical library 
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that could reproduce both empirical colour-temperature re- 
lations and, when combined with current isochrones, the ob- 
served colour-magnitude diagrams of star clusters. 

Thus, the present status of semi-empirical libraries is 
not yet fully satisfactory; it is still useful, however, to discuss 
these alternative libraries in connection with uncertainties 
in integrated magnitudes due to the spectral input. To this 
purpose we calculated an alternative set of SSPs based on 
the same isochrones as the G-SSPs (see §4) but adopting the 
semi-empirical library by Lejeune et al. (1997, 1998). We 
indicate this alternative set as L-SSPs; the corresponding 
solar magnitudes, derived for the solar model linked to the 
Lejeune library, are listed in Table 1. The difference between 
the G-SSPs and the L-SSPs is prominent in the V-I colour 
(Fig. 9) making the L-SSPs brighter in the I-band by up 
to 0.1 mag at intermediate and old ages (since the absolute 
V magnitude is no different between the two sets). This is 
purely an effect of the spectral library, and makes our model 
M/L ratios in the I band lower by a 5-10%, when the L-SSPs 
are adopted (Fig. 10). 

Finally, we checked whether a more refined treatment of 
the AGB phase would have a sizeable effect on the predicted 
M/L ratios with the Lejeune library, where red giants are 
more luminous in the I-band. The effect of a detailed inclu- 
sion of the TP- AGB phase, from the isochrones by Marigo 
& Girardi (2001), was found to be small in the I-band when 
the Girardi/Castelli library was adopted (§5.2 and Fig. 8). 
If the Marigo & Girardi isochrones are combined with the 
Lejeune library (set LM-SSPs), again the effects on the I- 
band M/L ratio is small with respect to the results from the 
L-SSPs (Fig. 10). In conclusion, with present models the I- 
band luminosity is more affected by the adopted spectral 
library than by the detailed modelling of the AGB phase. 

Since the situation with semi-empirical libraries is still 
somewhat dichotomic (Westera et al. 2002), in the remain- 
der of the paper we will adopt as our standard set the 
G-SSPs based on the theoretical spectral library by Gi- 
rardi/Castelli. However, from the discussion above and from 
Figs. 8 and 10, it is apparent that current alternative spec- 
tral libraries (and, marginally, a more refined modelling of 
the AGB phase) tend to increase the predicted I-luminosity. 
Therefore, our predictions for M/L ratios in the I band will 
be conservatively high, though not by large (at most 10%). 




10 9 
age [yr] 

Figure 9. Comparison between the G-SSPs and the L-SSPs 
(thick and thin lines, respectively) in V-I. The Lejeune library 
predicts redder V— I colours, hence more light in the I-band. The 
Kroupa IMF is adopted here for the sake of example, however the 
age evolution of V-I is insensitive to the IMF. 
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Figure 10. M/L ratios for our exponential disc models obtained 
with the L-SSPs (thin lines) as compared to the "standard" 
G-SSPs (thick lines). Dotted lines for the LM-SSPs. The Kroupa 
IMF is adopted here for the sake of example. 



6 CHEMO-PHOTOMETRIC MODELS 

More realistic photometric models for galactic discs should 
include star formation histories and corresponding chem- 
ical evolution consistently. From chemical evolution stud- 
ies of the Solar Neighbourhood, especially from the G- 
dwarf problem (Lynden-Bell 1975; Tinsley 1980; Pagel 1997 
and references therein), as well as from dynamical models 
of the formation of galactic discs (Larson 1976; Sommer- 
Larsen 1991; Burkert, Truran & Hensler 1992; Sommer- 
Larsen et al. 2003) it is well-known that galactic discs form 
gradually by slow accretion ("infall") of primordial or low- 
metallicity gas. Correspondingly, the SFR is expected first to 
increase, following the increase of available gas, reach a max- 
imum and later decrease due to gas consumption (Fig. 11). 

Moreover, the metallicity and colour gradients observed 
in the Galactic disc as well as external disc galaxies, in- 



dicate that the star formation and chemical enrichment 
history typically proceeds at a different (slower) pace at 
increasing galactocentric radius (inside-out scenario, Mat- 
teucci & Francois 1989; Chiappini, Matteucci & Gratton 
1997; Fig. 11; but see also Sommer-Larsen et al. 2003). A 
realistic SFH for a disc galaxy is therefore more complex 
than the monotonic, exponential form adopted in the sim- 
ple models of §5. 

As a base for our photometric models, we develop a 
series of chemical evolution models with infall, and with a 
multi-zone radial structure so that gradients of metallicity 
and SFH can be reproduced. The reader not keen on the 
details of model construction and calibration onto the ob- 
servable constraints, nor on the detailed discussion of the 
sets of models with each specific IMF, can directly find an 
overview of the results in §7.5. 

We make use of the chemical evolution code by Porti- 
nari, Chiosi & Bressan (1998); Portinari & Chiosi (1999). 
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Figure 11. Star Formation History for the disc of model SalpA 
(between R = 0.5 — 3 Hb, see §6.3). For comparison, the dotted 
line shows the SFH of a simple exponential model with the same 
6=0.79 and age 10 Gyr. The long-dashed and short-dashed lines 
illustrate the different pace of the SFH at different galactocen- 
tric distances, characteristic of inside-out models. For the sake of 
comparing the respective shape, all the SFR curves are normal- 
ized to the same integrated SFH of 1 Mq. 



The disc is assumed to have an exponential surface density 
profile: 



a(R) — a d e 



(3) 



where a(R) is the present-day total surface density distri- 
bution, Rd is the scale-length of the mass distribution and 
ad is the surface density at a distance R = Rd from the cen- 
ter. The disc is typically divided into 26 zones, 0.2 Rd wide, 
between R = and R = 5 Rd- 

In each zone infall is assumed to occur, as customary, 
at an exponentially declining rate: 



da 

— oc exp 
at \ t. 



t 



inf 

where the infall timescale may vary with radius, 
Tinf = Ti n f(R), typically increasing outward in the inside- 
out scenario. A similar behaviour of the infall process, expo- 
nentially declining rates with timescales increasing outward, 
is also found in some recent cosmological/hydro-dynamical 
simulations of the formation of galactic discs (Sommer- 
Larsen et al. 2003). 

We assume T = 13 Gyr for the age of the models, as sug- 
gested by the ACDM simulations of Sommer-Larsen et al., 
where the SF activity associated to the main galaxy typi- 
cally starts at a redshift of 4-5. This is larger than the age 
of 10 Gyr assumed in our exponential models in §5, but one 
should keep in mind that the latter represents the age of the 
bulk of local disc stars. In an infall scheme, the formation 
of the disc is described since the very early stages of gas 
accumulation; SF develops gradually, reaches its maximum 
after a few Gyrs and then declines, so that the the bulk of 
stars form in fact during the last ~10 Gyr (Fig. 11). 

As to the SF law, we adopt the form suggested by Do- 
pita & Ryder (1994): 

da, 1/3 5/3 
—jj- ct a 1 aj 

depending both on the gas surface density a g and on the 



total surface density a; a, is the surface density of stars. 
This is a Schmidt-like law (Schmidt 1959) with an efficiency 
decreasing at increasing radius through the dependence on 
a(R), as supported by both theoretical and empirical argu- 
ments (Dopita & Ryder 1994). The resulting SFHs are nec- 
essarily faster in the inner regions and slower in the outer 
ones where the surface density is lower; this SF law results 
in metallicity gradients comparable to the observed ones 
(Portinari & Chiosi 1999). With a convenient normalization 
(see Portinari & Chiosi 1999 for details) we express this law: 



d / j \ a 



1/3 (R,t) al /3 (R,t) 



a(Rd,T) 

where v is the star formation efficiency in units of \t~ x \. The 
code performs detailed calculations of the delayed gas and 
metal restitution from stars of various lifetimes, including 
supernovae of type la, and of the evolution of the abun- 
dances of the main chemical elements; we refer to Portinari 
et al. (1998) for all details on the chemical model and to 
Portinari & Chiosi (1999) for the implementation of the SF 
law and of the multi-zone scheme. With respect to the above 
mentioned papers, the present models are updated in the in- 
put yields of low and intermediate mass stars (Marigo 2001), 
which are consistent with the TP-AGB calculations included 
in the M-SSPs (§5.2). 

Our models are a two-parameter family, with parame- 
ters: 

Ti n f infall timescale (in Gyr; not to be confused with 
the SF timescale tsf of the simple models in §5) 
v SF efficiency (in Gyr" 1 ) 

This pair of parameters will be tuned, for each selected IMF, 
so as to reproduce two main constraints discussed below, the 
first selecting the relevant SFHs, the second setting the level 
of metal enrichment and the gas fraction. 



6.1 Constraints on the SFH 

From the SFH of each model we can define its 6-parameter 
V(T) 



6 = 



<ip> 



r = I J r[t),lt (I) 



where ^(t) is the instantaneous global SFR. We select the 
models with 6=0.8-1 as representative of Sbc-Sc discs, as 
discussed in §2; the corresponding typical intrinsic B-V 
colour is (B-V)o ~0.55 (Roberts & Haynes 1994). 



6.2 Constraints on chemical evolution 

A second constraint is the "typical chemical properties" of 
an Sbc-Sc disc. Non-barred spirals exhibit an oxygen gra- 
dient in the gas component that seems to be roughly inde- 
pendent of Hubble type, once radii are expressed in units of 
scalelength (in the B-band) rather than in physical distance 
units like kiloparsecs. The typical average value is: 



_d_ 



-0.2 dexh~ r 



with values ranging between -0.1 and -0.3 dex h B x (Garnett 
et al. 1997; van Zee et al. 1998; Prantzos & Boissier 2000). 
Oxygen can be taken as tracer of the global metal content, 
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as it is the most abundant metal and it generally contributes 
roughly half of the overall metallicity. In the following the 
term "metallicity" will be often used in place of "oxygen 
abundance" , for the sake of simplicity. 

The typical disc gas metallicity, usually defined as the 
central value extrapolated from the observed gradient or, 
more directly, the metallicity at a galactocentric distance 
R = fiB, shows a marked dependence on mass. The sam- 
ples by Garnett et al. and van Zee et al., as well as those by 
KTC94 and Roberts & Haynes (1994) considered above in 
the paper, mostly consist of galaxies with magnitude around 
Mb — —21 (when all rescaled to the same i/o=50). There- 
fore we will discuss mainly bright galaxies with sizes compa- 
rable to the Milky Way, while we do not investigate in this 
paper mass-dependent behaviours. The dependence of the 
typical metallicity on Hubble type seems to be a secondary 
relation, mostly a consequence of the mass-metallicity re- 
lation through the Hubble type-mass relation; there is in 
fact a tendency for earlier Hubble types to be more massive, 
though with a large overlap in mass between types. 

For Milky- Way sized galaxies, from the samples of Gar- 
nett et al. and van Zee et al. we derive a typical metallicity 
of Sbc-Sc discs of: 

Mi) 



12 + 1 



9.1 dex 



at R = Jib 



We will use this typical metallicity and metallicity gradient 
to constrain our chemical evolution models. 

Another significant constraint for chemical models is 
the gas fraction, that is the final ratio between the gas mass 
and the total baryonic mass (Tinsley 1980; Pagel 1997) 



fgas — 



M gas 



M gas + M* 



in obvious notation. Hence we need an estimate of the typ- 
ical gas fraction in Sbc-Sc discs. This in turn requires an 
estimate of the stellar mass M* , which can only be obtained 
from the observed luminosity and an assumed stellar M/L 
ratio (one cannot use dynamical mass estimates because of 
the unknown contribution of dark matter). However, since 
the M/L ratio is supposed to be a prediction of our models 
we should make no a priori assumptions about it, and in 
place of the usual gas mass fraction we consider as a con- 
straint the gas-to-luminosity fraction 

, _ M gas 

JL,gas — T 

Lib 

This quantity, besides being independent of any assumed 
M/L ratio, has the additional advantage of being inde- 
pendent of distance and hence of Ho. As the mass of gas 
(HI+H2) is derived from the observed flux in the 21 cm line 
and in the CO lines, which have the same dependence on the 
assumed distance as Lb, their ratio is distance independent. 

Roberts & Haynes (1994) give for Sbc-Sc galaxies a 
typical value of 



Lib 



Me 



Comparable values can be derived from the samples of de 
Mello et al. (2002a,b) and just slightly lower ones from that 
of KTC94. The masses of neutral hydrogen quoted here are 
the global masses for the whole HI disc. For the sake of chem- 
ical evolution, however, the significant constraint is the gas 



fraction within the area of active star formation, i.e. limited 
to within the optical radius. Some corrections to the above 
quantities should be applied to account for the fact that 
the HI disc generally extends farther than the optical disc. 
Considering the sample of spirals of KTC94, it is possible 
to correct from the total gas mass to the gas mass within 
the optical radius, using the "Roberts times" listed in their 
Table 5. We obtain typical corrections of 0.6 for field spi- 
rals and 0.8 for cluster spirals, whose outer gaseous discs 
are presumably stripped because of the dense environment, 
and hence less extended beyond the optical disc. Since the 
large sample of Roberts & Haynes, from the RC3 catalogue, 
contains both field and cluster spirals, we apply an aver- 
age correction of 0.7 to derive the gas masses within R 2 5- 
This agrees also with the corrections of 0.66-0.72 derived 
by Sommer-Larsen (1996). 

As to the molecular component, Young & Knezek 
(1989) give for Sbc-Sc galaxies 



M{H 2 ) 



= 0.16^°? 



Mq 



We find similar values in the sample of KTC94 and some- 
what lower values (~0.1) in the sample by of de Mello et al.^ 
Adding the two components, we derive: 

mi±m =0.3-0.37 

Lb 

and including the standard factor of 1.4 for the contribution 
of other elements (mostly helium) to the total gas mass, we 
finally derive 

f L ,gas ~ 0.42 - 0.52 

as a typical gas-to-luminosity fraction for Sbc-Sc galaxies. 

This value is broadly compatible with the analogous 
quantity in our Milky Way Disc. The masses of neutral 
and molecular hydrogen in our galaxy, within the opti- 
cal radius of the disc, are M(HI) = 4.3 x 10 9 M Q and 
M(H 2 ) = 1.3 - 1.7 x 10 9 M© (Binney & Merrifield 1998; 
Dame 1993). We can estimate the total B-band luminos- 
ity of the Milky Way as follows: the K-band luminosities of 
the Disc and Bulge are, respectively, Lk,t> = 4.9 x 10 10 Lq 
and L k ,b = 1-1 x 10 10 L (Binney & Merrifield 1998; Kent, 
Dame & Fazio 1991); considering (B-K)=3.51 as typical for 
a Sbc-Sc spiral (de Jong 1996b), we get for our Milky Way a 
global Lb ~ 1-8 x 10 10 Lq. Including the standard 1.4 fac- 
tor for the helium correction to the gas mass, we derive for 
the Milky Way f L , gaa ~ 0.45. 

Finally, as we model the SFH and the chemical evo- 
lution of galactic discs, we should adopt for the gas-to- 
luminosity fraction the luminosity of the disc only. So, we 
must correct for the contribution of bulges to the total B- 
luminosity. For Sbc-Sc galaxies, the average bulge contri- 



t Recently Boselli, Lequeux & Gavazzi (2002), applying a cal- 
ibrated conversion factor \ between CO emission and H 2 mass 
depending on metallicity, UV field and galaxy luminosity, derived 
a much lower typical M(H 2 ) / 'M(H I) <~ 0.15. They suggest that 
previous studies, adopting a constant conversion factor, overesti- 
mated the molecular mass in bright galaxies, while underestimat- 
ing it in faint galaxies, by a factor of 2-3. If so, our estimates of 
the gas-to-luminosity fractions above are upper limits, and the 
bulk of the gas would lie in the HI component. 
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bution is estimated to be Lb,b/Lb ~10-20% (Simien & de 
Vaucouleurs 1986; de Jong 1996a), which implies 



L,gas,T> 



0.45 - 0.55 



We therefore adopt 0.5±0.05 as our constraint for the chemi- 
cal enrichment level. However, this "gas fraction" constraint 
is much less firm than the observed metallicities, for it has 
been derived by assembling a number of independent obser- 
vational studies and applying various corrections. 



6.3 Calibration of the models 

For each of the six considered IMFs, we develop a grid of 
models with different timescales 7m/- For each model, the 
SF efficiency v is fixed so that the final gas metallicity is 12+ 
log(OfH) ~9.1 dex at T=13 Gyr and at a radius R = hs 
(see the discussion in §6.2 and Fig. 13). Among these models 
with specific {y,Ti n f) combinations, we select those whose 
global SFHs correspond to 6=0.8-1, that is to Sbc-Sc discs. 
This corresponds to selecting tendentially models with long 
infall timescales, as larger rj n /'s induce longer SFHs and 
hence larger b values. The parameters and characteristics of 
the selected models are listed in Table 2. 

The adopted SF law often generates metallicity gradi- 
ents in broad agreement with the observed ones, for it results 
in an inside-out SFH even when r OT / is constant with radius. 
Only when necessary, we include a radial increase of the in- 
fall timescale Tin/ (7?) reinforcing the inside-out behaviour 
(see model SalpC in Fig. 14). 

To illustrate our procedure, here we describe in detail 
the calibration of the models with the Salpeter IMF. The 
comparison to the results for the other "bottom-light" IMFs 
will be discussed in the next Section. 

Model SalpA adopts n n f = 5 Gyr, and the SF efficiency 
is fixed to be v = 0.47 Gyr -1 to get an oxygen abundance of 
9.1 dex at the reference radius R — hs (Fig. 13). We calcu- 
late the disc photometry by applying Eq. 2 to each individual 
annulus with its specific SFH ty(R,t) and metal enrichment 
history Z(R,t), to obtain the surface brightness profiles in 
Fig. 12. With an underlying exponential mass distribution, 
the light profile resulting from chemo-photometric models 
is not perfectly exponential, but close enough that a lumi- 
nosity scalelength can be defined at least within ~ 2 — 3 Rd 
(see also Boissier & Prantzos 1999 for similar results). In 
model SalpA, we derive a B-band scalelength h B ~ 0.95Rd 
within the inner ~ 3Rd- Fig. 12 also shows that the redder 
the band, the shorter the corresponding scalelength, as ex- 
pected and as observed also in the Milky Way (Boissier & 
Prantzos 1999 and references therein). 

Computing hs is necessary to define the limit of the 
optical disc, which is the region to be compared to ac- 
tual observations. A typical central (extrapolated) value of 
the surface brightness of discs is /iso=21.7 mag arcsec -2 
(Freeman 1970), although with a large scatter between 20 
and 22 mag arcsec -2 (for Hubble types up to Sc, de Jong 
1996a); the optical radius R25 corresponds to ^b=25, or 
R25 ~ 3 - Ah B from fj, B (R) = Mso + 1-086 R/h B . For the 
following, we define operationally the optical disc as the re- 
gion within R pt = 3 /is; we verified that taking 3.5 or 4 h b 
as the limit would not influence the global M/L ratios or gas 
fractions, since the outskirts of the disc are characterized by 
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Figure 12. Surface brightness in B, I and K-band for model 
SalpA, all normalized to /i = at _R = Rd- The thick line marks 
the B— band profile in the significant "disc region" between the 
bulge and the optical radius R op t (see text). Also shown is the 
underlying exponential mass profile cr(R), vertically displaced by 
-0.5 for clarity. 
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Figure 13. Metallicity gradient in the present-day gas phase 
for model SalpA, with radius in units of B-band scalelength h B . 
The significant "disc region" is marked as a thick line. The solid 
dot marks the observational calibration point 9.1 for metallicity 
at R = Kb- The dashed line is the observed average gradient 
-0.2 dex/h B , the dotted lines mark the limits of the broad range 
of observed gradients, between -0.1 and -0.3 dex/h B . 



a low surface density and brightness, and do not contribute 
significantly to the global values. 

We also exclude from our analysis the innermost 0.5 hs, 
where the bulge typically lies in real disc galaxies. 

Also the metallicity constraints derived from observa- 
tions are expressed in terms of the B-band scalelength hs, 
rather than the scalelength Rd of the underlying mass dis- 
tribution (Eq. 3). In Fig. 13 the oxygen abundance profile 
is plot as a function of R/h B , like actual observations; the 
model gradient is well within the typical observed range. 

Longer infall timescales tend to produce shal- 
lower metallicity gradients at the present time. With 
Tinf = 10 Gyr, the gradient is very shallow if n n f is constant 
over the disc; we obtain a gradient of about -0.2 dex//is 
only by adopting a radially increasing infall timescale (model 
SalpC, Fig. 14; Table 2). 
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Table 2. List of models with calibrated parameters (columns 3 to 5) and results (columns 6 to 10). (1) Model name. (2) Adopted IMF 
with corresponding mass limits. (3) Infall timescale r in f at R = hg, in Gyr. (4) SF efficiency v, in Gyr^ 1 . (5) Radial dependence of 
the infall timescale. (6) Scale-length of the B-band profile. (7) b-paramctcr of the SFH (Eq. 4). (8) M/L ratio in the I— band for the 
stellar component (stars+remnants). (9) Intrinsic (B-V) colour of the disc (including the central regions, i.e. between R = — 3 hg). 
(10) Gas-to-luminosity fraction. 
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Figure 14. Metallicity gradient for model SalpC, with radially 
increasing infall timescale (Table 2), compared to a model with a 
constant infall timescale r in f = 10 Gyr (dashed line). 

7 MODEL RESULTS 

In this Section we present the resulting M/L ratios and 
chemical properties of our full chemo-photometric models. 
First we discuss results for the Salpeter IMF, then for the 
other, "bottom-light" IMFs. All models are calibrated to 
reproduce the observed metallicity distributions and SFHs, 
as described in §6.3; the calibrated parameters and model 
results are all listed in Table 2. 

7.1 Salpeter models 

The I-band M/L ratios for the disc models SalpA-B-C, 
adopting the Salpeter IMF with the standard mass limits 
[0.1 — 100] Mq, are shown as solid circles in Fig. 15 as a 
function of the corresponding 6-parameter. Also shown as a 
continuous thin line is the prediction from the simple expo- 
nential models of §5, with solar metallicity. The full chemo- 
photometric models confirm quite well the high M/Li ~ 1.5 
predicted by the simple models for the Salpeter case. 

One should not infer, however, from the agreement with 
the simple models that photometric computations at con- 
stant solar metallicity (like those in §5) are sufficient. In fact, 
if we calculate the luminosity of the full chemical models 
SalpA-B-C using only SSPs with Z =Zq, we obtain higher 
M/L ratios (open circles in Fig. 15), mainly because of a 
slightly larger average age of stars in realistic models (with 
T = 13 Gyr) with respect to the simple exponential ones 
(with T = 10 Gyr). In full chemo-photometric models, this 
is compensated by the fact that an important fraction of the 
long-lived stars contributing to the I-band luminosity forms 
at metallicities below solar, lowering the M/L ratio with re- 
spect to calculations at fixed Z =Zq (solid circles vs. open 
circles; see also the M/L ratios at Z=0.008 in Fig. 4). So, 
the agreement with the simple exponential models with solar 
metallicity is the complex result of combining a more real- 
istic, radially dependent SFH, with metallicity-dependent 
photometric calculations. 

The crosses in Fig. 15 represent the result of calculating 
the photometry of models SalpA-B-C including only SSPs 
of metallicity Z <Zq, adopting SSPs of Z=Zq for larger 
metallicities. The difference with respect to the reference 
models is marginal, indicating that super-solar and solar 
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Figure 15. M/L ratio in the I— band for chemo-photometric mod- 
els with the Salpeter IMF. The shaded area marks the range 
M/Lr=0.7— 1 favoured by observations of Sbc-Sc discs (b = 0.8 — 
1). Solid circles: models SalpA-B-C adopting the standard IMF 
mass range [0.1—100] Mq. Open circles: models SalpA-B-C with 
photometry calculated at constant metallicity Z = Zq = 0.019. 
Crosses: models SalpA-B-C with photometry calculated including 
only SSPs with Z < Zq. Triangle: model SalpD with a smaller up- 
per mass limit M s = 70 Mq. Square: model SalpE with a smaller 
lower mass limit Mi = 0.09 Mq . Thin line: prediction from simple 
exponential models with Z = Zq. 




1 _ Salpeter 

0.6 0.7 0.8 0.9 l 

b 

Figure 16. M/L ratio in the B-band for chemo-photometric 
models with the Salpeter IMF. Symbols as in Fig. 15. 

metallicity stars contribute to the I-luminosity in a compa- 
rable way. However, neglecting super-solar metallicities in 
the photometry has noticeable effects when considering, for 
instance, the B-luminosity (Fig. 16). This is because the B- 
band is both more sensitive to metallicity than the I-band 
(Fig. 2) and more influenced by young stars, which have 
a metallicity that is higher on average and reaches super- 
solar values in the inner disc regions. Therefore, we stress 
the importance of covering the full metallicity range, from 
sub-solar to super-solar metallicities, when computing the 
photometric properties of disc galaxies. 

Fig. 17 shows the resulting gas-to-luminosity fraction 
Mg aa /LB for the chemo-photometric models, compared to 
the observational estimates for Sbc-Sc discs (§6.2). Models 
with the Salpeter IMF between [0.1-100] Mq (solid circles) 
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Figure 17. Gas-to-luminosity fraction for chemo— photometric 
models with the Salpeter IMF. Symbols as in Fig. 15. The hori- 
zontal lines mark the observational range of 0.5± 0.05 for the gas 
fraction in Sbc-Sc discs (6 = 0.8 — 1). 



tend to be too rich in gas with respect to observations. Con- 
sidering that for a given set of stellar yields, the final metal- 
licity increases with decreasing gas fraction (Tinsley 1980), 
our results indicate that with this IMF the global metal 
production (yield) is somewhat high: if we had imposed the 
models to reach the observed gas fractions, the final oxygen 
abundance in the gas would have been higher than observed. 
This finding is in line with known results for the Solar Neigh- 
bourhood: a Salpeter IMF extended up to 100 Mq produces 
too much oxygen to be compatible with the local evidence, 
while better agreement is found with an upper mass limit of 
40-50 M Q (Tsujimoto et al. 1997; Thomas, Greggio & Ben- 
der 1998; Gratton et al. 2000). 

The efficiency of metal enrichment from a stellar pop- 
ulation can be estimated by computing the so-called global 
(or net) yield (Tinsley 1980; Pagel 1997): 



yz 



pz(M)<f>(M) dM 



(5) 



M(T) 



where M(T) is the stellar mass of lifetime corresponding to 
the age T of the system (the galaxy), and pz(M) is the mass 
fraction of new metals ejected by a star of mass M; hence the 
integral expresses the amount of metals globally produced by 
a stellar generation over the age of the galaxy. The symbol 
a represents the fraction of mass locked up in "ever-living" 
low mass stars or remnants. In the Instantaneous Recycling 
Approximation, the net yield can be directly used to pre- 
dict the metallicity evolution (see Tinsley 1980; Pagel 1997 
for details); with more complex chemical networks taking 
into account finite stellar lifetimes, delayed gas restitution 
and the different timescales of release of different elements, 
as in our models, things are not as straightforward but the 
net yield per stellar generation still provides a useful in- 
sight on the efficiency of metal enrichment. In particular it 
underlines that the latter does not depend merely on the 
amount of metals produced per mass involved in star for- 
mation, but on the ratio between this and the mass that 
remains forever locked in stars. The yield can thus be modi- 
fied by altering the amount of metals produced (for instance 
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Figure 18. Upper panel: M/L ratio in the I-band for chemo- 
photometric models with the Kennicutt IMF; lower panel: cor- 
responding gas-to-luminosity fractions. The shaded area in the 
top panel and the dashed lines in the lower panel mark the obser- 
vational ranges for Sbc-Sc discs (b=0.8— 1). Solid circles: models 
KennA-B-C adopting the standard IMF mass range [0.1—100] Mq. 
Open circles: models KennD-E-F adopting the lowest possible 
mass end, with range [0-100] Mq. Solid triangles: models KennG- 
H-l-J with mass limits tuned to reproduce the observed gas frac- 
tions (Mi=0.05-0.1 Mq, M s =30-35 Mq). 



reducing/increasing the number of massive stars) and/or by 
changing the locked-up fraction. 

For the Salpeter models, a lower yield seems to be neces- 
sary since oxygen is somewhat over-produced. We calculate 
model SalpD decreasing the upper mass end M s of the IMF 
to reduce the metal production, and obtain a corresponding 
gas fraction in agreement with observations (solid triangle 
in Fig. 17 and 15). Alternatively, in model SalpE we decrease 
the lower mass end Mi and thereby increase the locked-up 
fraction, decrease the net metal yield and obtain again a gas 
fraction close to the observed values (solid square in Fig. 17). 

Notice however that, with the Salpeter IMF, a very 
restricted range of parameters (mass limits and infall 
timescales) can yield a Sbc/Sc-type SFH (b ~ 0.8) together 
with metallicities and gas fractions both close to the observa- 
tional values. What matters for the metallicity-gas fraction 
relation is in fact the yield, and models SalpD and SalpE have 
a similar oxygen yield (Table 3). However, when this condi- 
tion is met by decreasing Mi (model SalpE), the M/L ratio 
increases for we are shifting the IMF toward star masses that 
contribute no light (model SalpE, solid square in Fig. 15). 
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Figure 19. Upper panel: M/L ratio in the I— band for chemo- 
photometric models with the Larson IMF; lower panel: corre- 
sponding gas— to-luminosity fractions. The shaded area in the top 
panel and the dashed lines in the lower panel mark the observa- 
tional ranges for Sbc-Sc discs (6=0.8—1). Circles: models LarsA-B 
adopting the IMF mass range [0.01-100] Mg. Triangles: models 
LarsC-D-E with upper mass limit tuned to reproduce the observed 
gas fractions (M s =22-23 Mq). 



Figure 20. Same as Fig. 18 to 22, but for the Chabrier IMF. 
Circles: models ChabA-B-C adopting the full IMF mass range 
[0.01-100] Mq. Triangles: models ChabD-E-F with upper mass 
limit M s =32-33 Mq, tuned to reproduce the observed metallici- 
ties and gas fractions with high 6-parameters. 



M/L ratios are low, M/L/=0.7-0.8, still in agreement with 
the observed range. 



7.2 Kennicutt models 

The I-band M/L ratios of models KennA-B-C, adopting 
the Kennicutt IMF and standard limits [0.1 — 100] Mq, are 
shown as solid circles in Fig. 18 (upper panel); they lie close 
to the predictions from the simple exponential models of §5 
(thin line) and within the range favoured by observations. 

The lower panel shows that the corresponding gas-to- 
luminosity fractions M ga3 / Lb are higher than observed, by 
a factor of two or so. Hence, the net yield of this IMF is 
also too high (see the discussion in the previous section for 
the Salpeter models). Although the slope at the high-mass 
end (-1.5) is steeper than in the Salpeter IMF, so that one 
might expect less massive stars and less oxygen production, 
for the bottom-light Kennicutt IMF the smaller number of 
low-mass stars results in a low locked-up fraction and hence 
a high net yield overall (Table 3). The slope at the low-mass 
end is so shallow that even extending the IMF down to Mi=0 
would increase the locked-up fraction only marginally, and 
would not improve much the situation (models KennD-E-F, 
open circles in Fig. 18). To reduce the global yield, the oxy- 
gen production per stellar generation must also be reduced, 
by lowering the high-mass end. 

In models KennG-H-l-J (triangles in Fig. 18) we tune 
the mass limits of the IMF so that we obtain a SFH typical 
of Sbc/Sc discs (b — 0.8 — 1) together with a suitable gas- 
to-luminosity fraction M ga3 /LB ~ 0.5. The corresponding 



7.3 Larson and Chabrier models 

Models with the Larson or Chabrier IMF and mass limits 
[0.01-100] Mq (LarsA-B and ChabA-B-C in Table 2) have a 
low M/L ratio, in very good agreement with the correspond- 
ing simple models of §5; they are displayed as solid circles in 
Fig. 19 and 20, respectively. However, in both cases the cor- 
responding gas-to-luminosity fraction is exceedingly high, 
larger than the observed one by a factor of 2 for the Chabrier 
IMF and a factor of 4 for the Larson IMF. This again indi- 
cates that the global net yield is too high with these IMF 
(Table 3). This is no surprise for the Larson IMF which has 
a Salpeter slope at high masses — and the Salpeter IMF 
already had a somewhat high oxygen yield, §7.1 — and in 
addition, due to the exponential cut-off at low masses, its 
locked-up fraction is much smaller, leading to a very high 
net yield. The Chabrier IMF is steeper at high masses, yet 
the low locked-up fraction produces a high yield. 

Both IMFs have a sharp cut-off at low masses, so that 
even reducing the low mass end down to Mi = has no 
effect on the locked-up fractions, and the only way to reduce 
the net yield is to lower the upper mass end M s . Models 
LarsC-D-E with M s = 22 - 23 Mq, and models ChabD-E-F 
with M s =32-33 Mq reach the correct gas fraction, while the 
corresponding I-band M/L ratios are relatively unaffected 
by this tuning of the upper mass limits (triangles in Fig. 19 
and Fig. 20; see also Tsujimoto et al. 1997). 
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Figure 21. M/L ratio in the I-band for chemo-photomctric mod- 
els with the Kroupa IMF. The shaded area marks the range 
M/Lj=0.7— 1 favoured by observations of Sbc-Sc discs (b = 
0.8 — 1). Circles: models adopting the standard IMF mass range 
[0.1—100] Mq; they lie quite close to the predictions from the ex- 
ponential models from §5 (thin line). Triangles: models adopting 
a higher low mass end Mi=0.2— 0.4 to reproduced the observed 
metallicities at higher 6-values; the M/L ratio correspondingly 
decreases because of the lower locked— up fraction. Solid symbols: 
models KrouA-B-C-D-E listed in Table 2, with a gas-to— luminosity 
fraction M gas /LB ~ 0.5. Open symbols: other calculated models, 
with slightly different gas fraction. 

7.4 Kroupa and modified— Larson models 

The Kroupa IMF and the modified-Larson IMFs are char- 
acterized by a rather steep slope of -1.7 for M > 1 Mq. 
Combined with the locked-up fractions corresponding to the 
standard mass limits of [0.1-100] Mq, such steep slopes re- 
sult in low net yields compared to the IMFs previously con- 
sidered, especially for oxygen which is produced in massive 
stars (Table 3). — Actually, for the modified-Larson IMF 
the reference mass limits are [0.01-100] Mq, however due to 
the sharp cut-off below the peak mass ~0.25 Mq, changing 
the low-mass end of this IMF to Afj=0.1 (or 0) Mq has 
negligible effects on the results. 

As a consequence, with these IMFs the typical metallic- 
ities observed in disc galaxies can be obtained only for short 
infall timescales, and correspondingly low 6-parameters 
(models KrouA and LmodA-B-C in Table 2, displayed as solid 
circles in Fig. 21 and 22). Longer infall timescales tend in 
fact to dilute the metals produced and lower the final metal- 
licity, for a given net yield. 

To reach higher 6-values (b — 0.8 — 1), longer infall 
timescales are required; at the same time, the net yield must 
be increased to reproduce the typical oxygen abundance of 
9.1 dex at R = hs- This we can achieve by increasing the 
lower mass end Mi, which corresponds to decreasing the 
locked-up fraction and also the M/L ratio with respect to 
the standard case. 

As to the Kroupa case, models KrouB-C-D-E correspond 
to increasing values of the ^-parameter (solid triangles in 
Fig. 21). For simplicity, in Table 2 we list only Kroupa mod- 
els with a gas-to-luminosity fraction close to the observed 
value Mgas/LB ~ 0.5 (solid symbols in Fig. 21); other mod- 
els with slightly different gas fractions (open symbols) do 
not appreciably change the trends of the M/L ratios vs. 
6-parameter. A good Sbc/Sc model (b ~ 0.85) is KrouE, 
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Figure 22. Same as Figs. 18, 19 and 20, but for the modified 
Larson IMF. Circles: models LmodA-B-C adopting the full IMF 
mass range [0.01-100] Mq. Triangles: models LmodD-E-F with a 
lower mass limit of M;=0.2 Mq, tuned to reproduce the observed 
metallicities and gas fractions with high fe-parameters. 

but such high values of b, combined with the constraints on 
oxygen production, are reached only with an uncomfortably 
large Mi = 0.4 Mq. This limit, however, should not be taken 
at face value: a sharp cut off at low masses is just a simplistic 
way to reduce the locked-up fraction so as to increase the 
yield. Stars below ~0.8 Mq do not contribute metals nor 
substantial luminosity, and for chemical and photometric 
evolution their mass distribution is irrelevant, what matters 
is just their global mass. Any other distribution of masses at 
the low mass end, with the same locked-up fraction, would 
leave our results unchanged (see also Tsujimoto et al. 1997). 
Not surprisingly, with such a low locked-up fraction, the re- 
sulting M/L ratio is very low ( < 0.7, Fig. 21). 

For the modified-Larson IMF, the net yield is larger 
than in the Kroupa case, thanks to the exponential cut-off at 
low masses and the correspondingly lower locked-up fraction 
(for the same mass limits). In this case, a less extreme value 
of the low mass end M;=0.2 Mq suffices to reproduce the 
observed metallicities and gas-to-luminosity fractions with 
b > 0.8; the resulting typical M/L ratio is ~0.8. 

In conclusion, the Kroupa and modified-Larson IMFs, 
adopting a steep Scalo slope of -1.7 for M > 1 Mq, seem 
to meet the opposite problem than the IMFs previously dis- 
cussed: they are too inefficient in terms of oxygen production 
to account easily for the observed metallicities in spirals, if 
a SFH of the type of Sbc/Sc discs is required at the same 
time (b > 0.8). See however Kroupa (2001) about the un- 
certainties on the slope of the IMF at high masses: a value 
shallower than 1.7, possibly up to the Salpeter slope ~1.3, 
is by no means excluded from local observations. 
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7.5 An overview of the results 

For all the IMF considered, the full chemo-photometric 
models lie close to the predictions of the corresponding sim- 
ple models with the same mass limits — a complex result of 
combining realistic, infall-like and radially dependent SFHs 
with metallicity-dependent photometry. As a consequence, 
the conclusions of §5 for the M/L ratios are confirmed. 

The Salpeter IMF results in M/L ratios much higher 
than required to match observations; moreover, a Salpeter 
slope over the whole stellar mass range is currently not 
supported by observations (see §1). As to the bottom-light 
IMFs, broadly speaking we see two types of behaviour. 

In the first case the low locked-up fraction typical 
of bottom-light IMFs leads to a high global yield (Kenni- 
cutt, Larson and Chabrier models), too high to reproduce 
the typical gas fractions of Sbc/Sc spirals at the same time 
as their metallicities. These IMFs have yields comparable to 
what is required to explain the observed metal enrichment 
of hot gas in galaxy clusters (y ~2-3 Zq, Pagel 2002), and 
much higher than what is typically inferred for the Solar 
Neighbourhood (y < Zq) — compare the oxygen yields yo 
in Table 3 to the solar oxygen abundance Zo,q ~ 9 x 10~ 3 . 
A possible explanation is that disc galaxies, just like ellipti- 
cals in clusters, disperse a large part of the metals they pro- 
duce into the intergalactic medium, so that we do not detect 
these metals in present-day discs. Blow-out of material from 
present-day discs is observed in the form of fountains, but 
this gas does not seem to have the necessary energy to escape 
the galactic potential, so that it just falls back onto the disc 
on relatively short timescales (~ 10 s yrs, Bregman 1980; 
Fraternali et al. 2001; Heckman 2002). So, escape of large 
quantities of metals from present-day discs does not seem 
very realistic. However, in the early phases of galaxy forma- 
tion it is possible that dispersal of metals took place, when 
the potential well of the galaxy was weaker; large outflows 
are in fact observed in high-redshift Lyman-break galaxies, 
though these are usually identified with the progenitors of 
ellipticals (Pettini et al. 2000). Besides, there are arguments 
advocating large outflows of baryons from galaxies in general 
(Silk 2003). 

Alternatively, a large part of the metals produced may 
be locked-up in dust, which may correspondingly induce an 
observational underestimate of the true metallicity in spi- 
ral galaxies. Theoretical models and observational evidence 
suggest that in the interstellar medium of our Galaxy, 40% 
of the metals is locked-up in dust grains (Dwek 1998 and ref- 
erences therein). As a consequence, the oxygen abundance 
probed by atomic emission lines might be underestimated 
by almost a factor of 2. However, the oxygen abundance is 
typically measured in HII regions, where dust depletion is 
probably much lower than the average for the diffuse inter- 
stellar medium; besides, oxygen abundances and abundance 
gradients derived from the spectra of OB stars, which should 
be totally unaffected by dust depletion, agree very well with 
those from HII regions (Smartt & Rolleston 1997; Gummers- 
bach et al. 1998). So we consider oxygen abundances mea- 
sured in HII regions as a reliable estimate of the real metal- 
licity of a galaxy. Dust depletion does not seem to provide 
a viable solution to the apparent overproduction of metals. 

As an alternative solution to metal outflows, some fine- 
tuning of the upper mass limit of the IMF can be applied, 



corresponding to a tuning of the global yield, so that the 
chemical properties of discs can be matched without requir- 
ing considerable dispersal of metals. § However the resulting 
I-band M/L ratio is quite robust to such fine-tuning (see 
§7.6 and Tsujimoto et al. 1997), and for Sbc/Sc discs we 
consistently find for these IMFs M/Li < 1, with values be- 
tween 0.7 (Kennicutt IMF) and 0.9 (Chabrier IMF). 

In the second category, we find IMFs with a steep 
slope at the high-mass end (-1.7, Kroupa and modified- 
Larson IMFs) implying a low oxygen yield. In this case, 
to reproduce the observed metallicities and gas fractions 
combined with an Sbc/Sc-likc SFH (b >0.8) the global 
yield must be increased, via reducing the locked-up fraction 
with a higher minimum stellar mass M;. Correspondingly, 
since less mass is locked in very low-mass, low-luminosity 
stars, the M/L ratio also decreases, and for the Kroupa and 
modified-Larson models of Sbc/Sc galaxies we find again 
M/L/ -0.7-0.8. 

Our conclusions derived from the "chemical properties" 
of discs, i.e. the typical metallicities and gas fractions, are 
still subject to some uncertainties. The constraint of the gas- 
to-luminosity fraction is rather uncertain for it has been de- 
rived by assembling a number of independent studies and 
by applying various corrections (see §6.2); besides, it relies 
on the B luminosity, which is affected by uncertain dust ex- 
tinction. Furthermore, the global yield of a given IMF, and 
the corresponding fine-tuning of the mass limits necessary in 
chemical models, depends on the adopted stellar yields; here 
also some uncertainties exist, although the present study is 
mainly based on oxygen, which is the element with best 
estimated theoretical yields. However, the M/L ratio is gen- 
erally robust with respect to the fine-tuning of the IMF 
mass limits we applied to reproduce metallicities and gas 
fractions. For a metallicity distribution resembling the ob- 
served one (oxygen abundance of about 9.1 dex at R=hs, 
gradient around -0.2 dex/kpc) we get as a consistent re- 
sult that MjLi = 0.7 - 0.9 for Sbc/Sc disc models, with all 
the bottom-light IMFs we considered. Many other models 
were computed than those presented here, with shallower 
or steeper gradients (say, between -0.1 and -0.3 dex/kpc), 
or typical metallicities slightly different than the reference 
value of 9.1 (say between 9.0 and 9.2), and the conclusions 
on the M/L ratio in the I band for late-type discs are quite 
robust to such variations in the model constraints. 



7.6 The driving parameters for the M/L ratio 

We have demonstrated that various "bottom-light" IMFs 
suggested in recent literature result in low M/L ratios 
M/Li < 1 for Sbc/Sc discs. It is interesting to discuss what 
is the common feature, among these different IMFs, that 
determines this result. 

§ The same reduction of the net yield is obtained, if above a cer- 
tain stellar mass there is fallback of metals onto a central black 
hole after the supernova explosion (Colpi, Shapiro & Wasser- 
man 1996; Zampicri, Shapiro & Colpi 1998; MacFadyen, Woosley 
& Heger 2001; Podsiadlowski et al. 2002). In the presence of 
fallback, the actual amount of metals enriching the interstellar 
medium is less than what is formally predicted by supernova mod- 
els. See also Tsujimoto et al. (1997). 
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Figure 23. 3-D plot relating the M/Li ratio to the b-parameter and to the IMF "driving quantity" C[i-2]/ Cf ( see text). 
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Figure 24. Relation between M/Lj and the "driving quantity" 
f[l — 2]/a (see text), for models with b >= 0.7. Solid symbols 
for the "calibrated" models in Table 2, open symbols for other 
models. A given IMF-symbol can correspond to different values 
of the abscissa because the models explore different mass limits 
(cf. Table 2). 



We analyzed in our models the correlations between the 
predicted M/L ratio and various quantities characterizing 
the IMF, e.g. the mass fraction in stars around 1 Mq, or 
around 10 Mq, or the mass fraction in "ever-living" stars 
below 1 Mq, etc. The best correlation is found with the 
quantity C[i-2]/ a , where C[i— 2] is the mass fraction that a 
given IMF "assigns" to stars between 1-2 Mq and a is the 
mass fraction that remains locked in ever-living stars and 
remnants. Fig. 23 shows a 3-D plot of the resulting M/Li as 
a function of the 6-parameter and of the above mentioned 
quantity £[ 1 _ 2 ]/a characterizing the IMF. The plot shows 



the "calibrated" models in Table 2, as well as many other 
models with a variety of IMF mass limits, metallicities, gas- 
to-luminosity fractions etc., run during the calibration pro- 
cess. Besides the expected dependence on the fc-parameter 
amply discussed in this paper, the typical M/Li evidently 
decreases at increasing ^ 1 _ 2 \/a. Fig. 24 underlines the latter 
dependence by showing a "slice" of the 3-D plot of Fig. 23, 
corresponding to large 6-parameters b > 0.7 and projected 
on the ( plane. These correlations indicate that, 

in general, the M/Li ratio is lower when the IMF favours 
stars in the [1-2] Mq range versus stars of lower masses and 
remnants from more massive stars. 

This "driving quantity" £n_2]/a is relevant for the M/L 
ratio in red optical bands like the I band. At different wave- 
lengths, the luminosity is dominated by stars of different 
masses and the corresponding IMF — M/L ratio relation is 
driven by a different quantity. For instance, a similar analysis 
of the M/L ratio in the B band shows that this is influenced 
mostly by the mass fraction in the range [1-9] Mq, as ex- 
pected from the fact that B-band light is sensitive to the 
contribution of more massive stars. 

For the IMFs and the corresponding mass limits in the 
calibrated models of Table 2, we list the characteristic quan- 
tities C[i-2]> a etc. in Table 3. We list also the global oxygen 
yield yo (see Eq. 5) and the related IMF quantity £9 /a, 
where £g is the mass fraction in stars more massive than 
9 Mq — since it is these stars that contribute the bulk of 
the oxygen production. Actually, stellar yields depend on 
metallicity and this effect is included in our chemical mod- 
els (Portinari et al. 1998); the global yields and locked-up 
fractions listed for illustration in Table 3 are computed for 
solar metallicity. 

Traditionally, in chemical models the stellar IMF has 
been characterized by the quantity fi , that is the mass frac- 
tion in stars more massive than 1 Mq. This corresponds to 
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Table 3. Characteristic quantities for the IMFs adopted in the 
models of Table 2 (see text), a is the locked— up fraction, C[i-2] i s 
the mass fraction in the [1-2] Mq range, £g is the mass fraction 
in stars more massive than 9 Mq. yo is the global oxygen yield 
of the IMF. 



dividing the IMF into two mass ranges, the "ever-living" 
low-mass stars, locking up mass out of further chemical pro- 
cessing forever, and stars with a lifetime shorter than a Hub- 
ble time which recycle gas and metals; this description is 
useful and meaningful within the simplified Instantaneous 
Recycling Approximation for chemical evolution. A similar 
approach is applied also when deriving the IMF from the 
present-day mass function (PDMF), grossly separating ever- 
living stars, whose PDMF is the integrated result of the IMF 
over the history of the galaxy, from stars with finite stellar 
lifetimes assumed then to evolve "instantly" , whose PDMF 
is a direct signature of the IMF (Tinsley 1980; Scalo 1986). 

In the light of our results for the M/L ratio, we sug- 
gest a gross description of the IMF based on three mass 
ranges. (1) Ever-living stars (M < 0.8-1 Mq) contribute to 
the locked-up fraction; their contribution to luminosity be- 
comes important only for very old systems, where they are 
not much smaller than the typical "turn-off mass" . (2) Inter- 
mediate living stars (M ~l-9 Mq) contribute an important 
amount of recycled gas (for a Salpeter IMF, ~2/3 of the total 
returned gas fraction) and dominate the luminosity; depend- 
ing on the observational band considered, the "important" 
sub-range for the luminosity may change: e.g. [1-2] Mq for 
the I-band, [1-9] Mq for the B-band. (3) Short-lived stars 
(M >9 Mq) contribute part of the recycled gas, part of the 
locked-up fraction (through the remnants) and most of the 
metals (oxygen in particular); their luminosity contribution 
is dominant only for blue and UV bands, and/or for very 
young systems or systems dominated by a recent burst. We 
prefer here the terminology "ever, intermediate and short 



living" stars rather than "low, intermediate and high mass" 
stars, to distinguish this classification from the mass ranges 
classified on the base of stellar evolution properties (Iben 
1991; Chiosi, Bertelli & Bressan 1992). 

Broadly speaking, it is the relative proportion of the 
mass fractions in these three ranges that determines the ex- 
pected chemo-photometric properties of the system. Since 
for the I-band, the mass range dominating the luminosity 
differs from the mass range responsible for the oxygen pro- 
duction, it is easy to understand how we could to tune the 
yields through the IMF upper mass limits and reproduce the 
observed gas-to-luminosity fraction, with negligible effects 
on the M/L/ ratio (§7.5; see also Tsujimoto et al. 1997). 
However, this dichotomy may not hold when one consider 
other bands (for instance blue or UV bands more sensitive 
to the contribution of massive stars) or other chemical ele- 
ments (for instance helium, carbon, nitrogen and iron, pro- 
duced in non-negligible amounts by "intermediate living" 
stars — including SN la) . On the contrary, one can envisage 
a connection between the M/L ratio in certain bands and 
the production of certain elements. So for a full account of 
the chemo-photometric evolution of galaxies, it is important 
to maintain consistency between the computed luminosity 
and chemical outputs (Marigo & Girardi 2001). 

7.7 The total baryonic M/L ratio 

In this section we discuss the global baryonic M/L ratio, 
including the mass in gas, for those models of Sbc/Sc discs 
(large b values) that also comply with the constraint on the 
gas-to-luminosity ratio (§6.2) — and hence have the gas 
amount indicated by observations. 

With the Salpeter IMF, the typical global baryonic M/L 
ratio in Sbc/Sc discs is > 2 Mq/Lq (Table 4, column 5), 
and the mass in gas is about 25% of the total baryonic mass 
(column 6). With the bottom light IMFs, the gas contribu- 
tion is more important since less mass is stored in the stellar 
component, for the same luminosity. With the adopted gas- 
to-luminosity ratio, the gas mass typically accounts for 40% 
of the total baryonic mass in Sbc/Sc discs, and the global 
baryonic M/L ratio is 1.2-1.5 Mq/Lq. 

A gas fraction of 40% may seem exceedingly large, but 
here we are dealing with late-type spirals which are noto- 
riously gas-rich (Roberts & Haynes 1994; Sommer-Larsen 
1996). Besides, the gas fraction is a strong function of ra- 
dius: in the inner parts, say within R < lhs, our Sbc/Sc 
models typically have a gas fraction less than 20%, increas- 
ing at higher radius. So, while the gas mass is relevant in 
computing the global baryonic mass in the disc, in the inner 
regions the baryonic content is definitely dominated by the 
stellar component. 

Besides, the face-value gas fractions we find depend 
on the adopted gas-to-luminosity ratio, suggested from ob- 
servations yet rather uncertain (§6.2). This section is just 
meant to draw attention on the fact that, with bottom-light 
IMFs, the gas mass is likely to contribute significantly to the 
total baryonic mass in late-type discs, and it should not be 
neglected when estimating the global baryonic content. This 
may be crucial when discussing the physical nature of the 
Tully-Fisher relation, since over a wide enough dynamical 
range the rotation velocity appears to correlate better with 
the global baryonic (stars+gas) mass, rather than with the 
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(1) 

Model 


(2) 

IMF 


(3) 

b 


(4) 

M, 
Li 


(5) 

Mbar 

Li 


(6) 

M gaB 

M bar 


SalpD 
SalpE 


Salpcter [0.1-70] 
Salpetcr [0.09-100] 


0.78 
0.79 


1.58 
1.67 


2.09 
2.23 


0.24 
0.25 


KrouD 
KrouE 


Kroupa [0.35-100] 
Kroupa [0.4-100] 


0.73 
0.84 


0.73 
0.66 


1.12 
1.06 


0.35 
0.38 


KennG 
KennH 
KennI 
KennJ 


Kennicutt [0.1-30] 
Kennicutt [0.1-35] 
Kennicutt [0.05-35] 
Kennicutt [0.05-35] 


0.78 
0.97 
0.81 
0.90 


0.74 
0.69 
0.79 
0.76 


1.22 
1.21 
1.35 
1.23 


0.40 
0.43 
0.42 
0.39 


LarsC 
LarsD 
LarsE 


Larson [0.01-22] 
Larson [0.01-22] 
Larson [0.01-23] 


0.81 
0.83 
0.95 


0.78 
0.79 
0.74 


1.37 
1.28 
1.30 


0.43 
0.39 
0.43 


LmodD 
LmodE 
LmodF 


mod. Larson [0.2-100] 
mod. Larson [0.2-100] 
mod. Larson [0.2-100] 


0.86 
0.90 
0.93 


0.82 
0.80 
0.80 


1.36 
1.30 
1.26 


0.40 
0.39 
0.36 


ChabD 
ChabE 
ChabF 


Chabricr [0.01-32] 
Chabrier [0.01-32] 
Chabricr [0.01-33] 


0.83 
0.87 
0.92 


0.94 
0.91 
0.90 


1.53 
1.45 
1.42 


0.38 
0.37 
0.36 



Table 4. Total baryonic M/L ratio of the model Sbc/Sc discs 
from Table 2, including the mass in gas. (1) Model name. (2) 
Adopted IMF with corresponding mass limits. (3) 6-parameter 
of the SFH (Eq. 4). (4) M/L ratio in the I-band of the stellar 
component (comprehensive of living stars and stellar remnants). 
(5) Total baryonic M/L ratio including the gas mass: Mi, ar = 
M» + M gas . (6) Gas mass fraction in the baryonic disc. 



stellar mass only (McGaugh et al. 2000). When discussing 
the baryonic TF relation, the M/L ratio ascribed to the 
stars and, correspondingly, the relative contribution of the 
gas mass have important consequences for the derived slope 
(McGaugh 2001, 2003). Lower-than-Salpeter M/L ratios, as 
favoured in the present paper, imply a somewhat shallower 
TF relation (Bell & de Jong 2001; de Jong & Bell 2001). 



7.8 An indirect comparison to the Solar 
Neighbourhood 

Tsujimoto et al. (1997, hereinafter T97) performed an anal- 
ogous study to our present one, combining chemical evo- 
lution and M/L ratio constraints for the Solar Neighbour- 
hood to determine a three-parameters IMF. Although their 
study was similar in spirit to ours, it is not trivial to com- 
pare results, both because of substantial differences in the 
models (in stellar yields, SN la rate formalism, SF law and 
photometric computations) and because of the different ap- 
proach (detailed chemical patterns and star counts in the 
Solar Neighbourhood, versus global average properties of 
external Sbc/Sc disc galaxies). Still, a comparison can be 
interesting and we shall broadly outline it. 

T97 considered single-slope IMFs with varying mass 
limits so that their set of IMF parameters is {Mi, M s , x) 
where x is the power-law slope. A direct comparison with 
our results is possible for the "power-law IMFs" considered 
in this paper (Salpeter, Kroupa and Kennicutt): although 
the Kroupa and Kennicutt IMFs are not single-slope, the 
slope changes only below 1 M Q where what matters is ef- 
fectively the global mass in low-mass stars, rather than the 



detailed mass distribution (T97 and §7.4). For our "best 
models" of Sbc/Sc discs listed in Table 4, with the right 
gas-to-luminosity ratio (i.e. level of metal enrichment), we 
can compute the "equivalent" value of Mi that, for a single- 
slope IMF, would result in the same mass fraction below 
1 M e {Mi ,eq in Table 5). We also compare the mass frac- 
tions within the three significative mass ranges discussed 
in §7.6: "ever-living" {M < 1 Mq), intermediate-living 
{M = 1 — 9 Mq) and short-lived (M > 1 M ) stars. In Ta- 
ble 5 we list the various IMF parameters as constrained by 
T97 and in this paper. For brevity, we consider the highest 
values of M s and of Mi allowed in Fig. 3 of T97; considering 
their lower values somewhat increases the discrepancy with 
our results (see below), but does not change the outcome of 
this qualitative comparison. 

For each of the IMFs in Table 5, our results favour 
smaller low-mass fractions £<i, and higher fractions both 
in the intermediate [1-9] Mq and in the high > 9 Mq mass 
range. This implies both a higher metal enrichment (roughly 
driven by C>9/ a or C>s>/C<i) an d a lower M/L ratio (roughly 
driven by C[i-9]/C<i) m our models. The higher level of 
metal enrichment might be related to the fact that the gas- 
to-luminosity ratio of external spirals is certainly a more 
uncertain and loose constraint on chemical evolution than 
the detailed abundance patterns in the Solar Neighbourhood 
considered by T97. For instance, a very recent compilation 
(Bettoni, Galletta & Garcfa-Burillo 2003) seems to favour 
lower gas fractions {Mhi/Lb ~ 0.2 and Mh 2 /Lb ~ 0.1) 
than we adopted in §6.2, which would correspond to a lower 
efficiency of metal enrichment for external spirals. In the 
context of our models, this could be obtained via a different 
tuning of the upper mass limit M 3 of the IMF to reduce the 
characteristic yield, but as discussed in §7.5 and §7.6 the 
resulting M/L ratios would hardly change. 

As to the different M/L ratio between T97 and this pa- 
per, it is difficult to understand whether the discrepancy is 
due to model differences or to the different approach (So- 
lar Neighbourhood vs. external discs). The difference is not 
simply accounted for by different SFHs, since both the Solar 
Neighbourhood in T97 and the Sbc/Sc discs here have been 
modelled with the respective suitable SFHs. Possibly, while 
chemical evolution is much better constrained in the Solar 
Neighbourhood, the estimate of the local luminosity from 
star counts (T97) is more cumbersome than the study of the 
global luminosity in external galaxies. Besides, dark matter 
may have a minor, but non-zero contribution to the local dy- 
namical M/L ratio, easing the M/L constraints adopted by 
T97. Alternatively, the discrepancy might be resolved when, 
for the M/L ratio of external galaxies, allowance is made for 
the dispersion and the uncertainty in the zero-point of the 
TF relation (Appendix A). 

The issue deserves further investigation, but this is be- 
yond the purpose of this paper. Reasonably, one should first 
assess whether there is an effective discrepancy in the obser- 
vationally deduced M/L ratio and IMF, free from modelling 
issues. The two environments, the Solar Neighbourhood and 
outer disc galaxies, should be analyzed with the same models 
so that a significant comparison can be made. We postpone 
to future work the study of the M/L ratio and IMF in the 
Solar Neighbourhood with our chemo-photometric models. 
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Salpctcr 


1.35 


0.05 


50 


0.71 


0.21 


0.08 


0.09-0.1 


70-100 


0.62 


0.26 


0.12 


Kcnnicutt 


1.5 


0.14 


55 


0.66 


0.26 


0.08 


0.35-0.38 


30-35 


0.44 


0.45 


0.11 




1.6 


0.17 


58 


0.68 


0.26 


0.06 












Kroupa(*) 


1.7 


0.23 


(51 


0.65 


0.29 


0.06 


0.44-0.47 


100 


0.44 


0.46 


0.10 



Table 5. Comparison between IMFs as constrained by T97 and in the present paper (PST). Mi and M s are the mass lim- 
its of the IMF of slope x\ C<1> C[i— 9] an d C>9 are t ne corresponding mass fractions in the three representative mass ranges. 
M j e9 (PST) is the "single-slope equivalent" lower mass limit for our IMFs (see text). 

(*) Slopes x > 1.6 are formally ruled out by T97, because for steep slopes they find that no value of Mi can fulfill both the chem- 
ical and the M/L ratio constraints. Still, for qualitative comparison with our Kroupa models we consider their case x = 1.7 with 
Mi = 0.23, an intermediate value between those allowed by chemical and M/L ratio arguments (Fig. 3 in T97). 



8 THE BULGE CONTRIBUTION 

In this section we attempt to give a rough estimate of the 
contribution of the bulge to the global M/L ratio of disc 
galaxies. Bulges are tendentially older and dimmer than 
discs, although in some cases a continuity of properties 
(colours etc.) is seen between the bulge and the inner re- 
gions of the disc, especially for late-type spirals (Balcells 
2003 and references therein). A lower limit to the global 
M/L ratio is hence given by considering the central regions 
as the extrapolation of the disc properties to the center, 
which in included in our models. As described in §6.3, from 
the analysis of our disc models we had excluded the inner 
regions (R < 0.5 /is); including them yields the M/L ratios 
in parenthesis in column 4 of Table 6. We limit the list to 
models with large b parameters, as representative of Sbc/Sc 
spirals. The inclusion of the central disc regions induces a 
minor increase (at most 10%) in the global M/L ratio. In 
particular, for all the "bottom-light" IMFs considered, the 
M/L ratio remains below 1 for the Sbc/Sc models. 

An upper limit to the bulge M/L ratio is given by 
the characteristic M/L ratio of an old SSP, neglecting any 
possible tails of later star formation and younger stellar 
components — expected among others because of gas re- 
cycling from the first population, and detected for instance 
in the bulge of the Milky Way (Ng et al. 1995; van Loon 
et al. 2003). We consider SSPs of age 10 Gyr, as character- 
istic of the bulges of early-type spirals (Peletier et al. 1999); 
this is likely to overestimate the M/L ratio of the bulges 
of Sbc/Sc galaxies, which tend to be younger (Peletier 
et al. 1999; Carollo, Ferguson & Wyse 1999). We consider 
SSPs of solar metallicity for simplicity. In the Bulge of our 
galaxy, the bulk of stars is between half-solar and solar 
metallicity (Mc William & Rich 1994; Zoccali et al. 2003) 
with a super-metal rich tail (Bertelli et al. 1995); varying 
the metallicity by a factor of 2 would change the bulge M/L 
ratio in the I-band by 15% or so, negligible for the sake of 
the rough estimates we are giving here. Column 5 of Ta- 
ble 6 lists the I-band M/L ratios of the stellar component 
(i.e. stars+remnants, excluding the re-ejected gas) of 10 Gyr 
old, solar SSPs with the same IMF as in the disc models. 

Typically, the bulge-to-disc luminosity ratio in the I— 
band, (B/T>)j, is around 25% for Sbc/Sc spirals, although 
with a large scatter (Sommer-Larsen et al. 2003). Assuming 
(i3/Z?) j=0.25 and the respective M/L ratios for disc and 
bulge, results in mass ratios between 50 and 70% (column 6). 
This illustrates that the luminosity ratio, even in the red 
bands, is not a direct indicator of the mass ratio of the two 
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Figure 25. Estimated global M/L ratio of disc galaxies, including 
the contribution of the bulge. Dots represent the lower limit es- 
timate (treating the bulge as an extension of the disc properties 
to the centre), bars indicate the upper limit estimate (treating 
the bulge as a 10 Gyr old SSP contributing 25% of the total I— 
band luminosity). The shaded area marks the range M/L/=0.7— 1 
favoured by observations of Sbc— Sc spirals. 



components: allowance must be made for the different M/L 
ratio between an old bulge and the disc. 

With the adopted (B/2?) / =25% and the respective M/L 
ratios, we estimate the global galactic M/L ratios in col- 
umn 7 of Table 6. For the bottom-light IMFs, MjLi (tot) = 
1 — 1.2 (see also Fig. 25). These values are strict upper lim- 
its, for two reasons. Firstly, as mentioned above, bulges in 
Sbc/Sc spirals seem to be younger and/or host a younger 
component than the assumed 10 Gyr, so their character- 
istic M/L ratio is expected to be lower than in column 5. 
Secondly, we are likely overestimating here the luminosity 
contribution from the bulge: the typical (B/2?), ~0.25 in 
Sommer-Larsen et al. (2003) are in fact derived from 2D 
light decompositions assuming a r~ profile for the bulge 
(Byun 1992). Nowadays, exponential light profiles or shal- 
lower Sersic profiles are favoured for bulges, especially in late 
type spirals (de Jong 1996a; Carollo et al. 1999, 2001; Bal- 
cells 2003); with an exponential profile, the decomposition 
yields significantly lower B/T> ratios (de Jong 1996a) and the 
bulge contribution to the total M/L ratio is correspondingly 
less. Notice also that the global (B-V)o colours in Table 2, 
obtained for the global disc (extended to the center) are in 
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(i) 

Model 



(2) 

IMF 



(3) 

b 



(4 

M» 
Li 



V) 



(5) 



M, 
Li 



(B) 



(6) 
M, (B) 
M» (X>) 



(7) 



(8) 



(9) 



(tot) 



(D) 



M, 



(B) 



(10) 



SalpA 


Salpcter [0.1-100] 


0.79 


1.54 


(1.60) 


3.77 


0.61 


1.99 


0.67 


1.31 


0.31 


SalpB 


Salpeter [0.1-100] 


0.86 


1.52 


(1.57) 


3.77 


0.62 


1.97 


0.65 


1.31 


0.31 


SalpC 


Salpeter [0.1-100] 


0.90 


1.51 


(1.62) 


3.77 


0.62 


1.96 


0.64 


1.31 


0.30 


SalpD 


Salpeter [0.1-70] 


0.78 


1.58 


(1.68) 


3.76 


0.59 


2.02 


0.66 


1.31 


0.30 


SalpE 


Salpeter [0.09-100] 


0.79 


1.67 


(1.71) 


3.99 


0.60 


2.13 


0.70 


1.38 


0.30 


KrouD 


Kroupa [0.35-100] 


0.73 


0.73 


(0.77) 


1.53 


0.52 


0.89 


0.31 


0.53 


0.31 


KrouE 


Kroupa [0.4-100] 


0.84 


0.66 


(0.68) 


1.41 


0.53 


0.81 


0.28 


0.49 


0.31 


KennA 


Kennicutt [0.1-100] 


0.84 


0.75 


(0.77) 


1.80 


0.60 


0.96 


0.33 


0.61 


0.32 


KennB 


Kennicutt [0.1-100] 


0.93 


0.73 


(0.74) 


1.80 


0.62 


0.94 


0.32 


0.61 


0.32 


KennC 


Kennicutt [0.1-100] 


1.02 


0.71 


(0.72) 


1.80 


0.63 


0.93 


0.31 


0.61 


0.32 


KennD 


Kennicutt [0-100] 


0.82 


0.92 


(0.95) 


2.23 


0.61 


1.18 


0.35 


0.76 


0.28 


KennE 


Kennicutt [0-100] 


0.89 


0.90 


(0.93) 


2.23 


0.62 


1.17 


0.34 


0.76 


0.28 


Kenn F 


Kennicutt [0-100] 


1.00 


0.87 


(0.89) 


2.23 


0.64 


1.14 


0.33 


0.76 


0.29 


KennG 


Kennicutt [0.1-30] 


0.78 


0.74 


(0.76) 


1.78 


0.60 


0.95 


0.31 


0.61 


0.31 


KennH 


Kennicutt [0.1-35] 


0.97 


0.69 


(0.71) 


1.79 


0.65 


0.91 


0.30 


0.61 


0.32 


Kennl 


Kennicutt [0.05-35] 


0.81 


0.79 


(0.81) 


1.94 


0.61 


1.02 


0.34 


0.66 


0.31 


KennJ 


Kennicutt [0.05-35] 


0.90 


0.76 


(0.80) 


1.94 


0.64 


1.00 


0.32 


0.66 


0.31 


LarsA 


Larson [0.01-100] 


0.82 


0.84 


(0.86) 


2.16 


0.64 


1.10 


0.37 


0.74 


0.32 


LarsB 


Larson [0.01-100] 


0.94 


0.78 


(0.80) 


2.16 


0.69 


1.06 


0.34 


0.74 


0.32 


LarsC 


Larson [0.01-22] 


0.81 


0.78 


(0.84) 


2.11 


0.68 


1.05 


0.33 


0.72 


0.31 


LarsD 


Larson [0.01-22] 


0.83 


0.79 


(0.84) 


2.11 


0.67 


1.05 


0.33 


0.72 


0.31 


LarsE 


Larson [0.01-23] 


0.95 


0.74 


(0.79) 


2.11 


0.71 


1.01 


0.31 


0.72 


0.31 


LmodB 


mod. Larson [0.01-100] 


0.74 


0.93 


(0.99) 


2.18 


0.59 


1.18 


0.40 


0.75 


0.31 


LmodC 


mod. Larson [0.01-100] 


0.77 


0.95 


(1.01) 


2.18 


0.57 


1.20 


0.40 


0.75 


0.30 


LmodD 


mod. Larson [0.2-100] 


0.86 


0.82 


(0.85) 


1.96 


0.60 


1.05 


0.35 


0.75 


0.28 


LmodE 


mod. Larson [0.2-100] 


0.90 


0.80 


(0.84) 


1.96 


0.61 


1.03 


0.34 


0.68 


0.31 


LmodF 


mod. Larson [0.2-100] 


0.93 


0.80 


(0.84) 


1.96 


0.61 


1.03 


0.34 


0.68 


0.31 


ChabA 


Chabricr [0.01-100] 


0.78 


0.95 


(0.97) 


2.54 


0.67 


1.27 


0.41 


0.87 


0.32 


ChabB 


Chabrier [0.01-100] 


0.88 


0.91 


(0.93) 


2.54 


0.70 


1.24 


0.39 


0.87 


0.32 


ChabC 


Chabrier [0.01-100] 


0.95 


0.88 


(0.90) 


2.54 


0.72 


1.21 


0.38 


0.87 


0.31 


ChabD 


Chabricr [0.01-32] 


0.83 


0.94 


(0.97) 


2.53 


0.67 


1.26 


0.40 


0.87 


0.31 


ChabE 


Chabricr [0.01-32] 


0.87 


0.91 


(0.97) 


2.53 


0.70 


1.23 


0.38 


0.87 


0.31 


ChabF 


Chabrier [0.01-33] 


0.92 


0.90 


(0.96) 


2.53 


0.70 


1.23 


0.38 


0.87 


0.31 



Table 6. Estimated bulge contribution, total (M/L)j ratios and K-band B/T> ratios for the Sbc/Sc models of Table 2. (1) Model 
name. (2) Adopted IMF with corresponding mass limits. (3) b— parameter of the SFH (Eq. 4). (4) M/L ratio in the I-band of the stellar 
component in the disc; in parenthesis the M/L ratio when considering the disc extended to the galactic center. (5) Stellar M/L/ ratio of 
a SSP of solar metallicity and age 10 Gyr, representing the bulge, with the IMF indicated in column 2. (6) Actual mass ratio between 
the disc and a bulge with I-band luminosity 25% of that of the disc. (7) Global galactic M/L/ ratio assuming a 10 Gyr old bulge with 
(B/I')j=25%. (8) M/L ratio in the K— band of the stellar component in our disc models. (9) Stellar M/Lx ratio of a SSP of solar 
metallicity and age 10 Gyr, representing the bulge, with the IMF indicated in column 2. (10) B/V luminosity ratio in the K-band 
corresponding to the mass ratio in column 6, or equivalcntly to the adopted (Z3/X>) / =25% 



agreement with the observed (B-V)~0.55 of Sbc/Sc 
spirals (Roberts & Haynes 1994); hence the contribution of 
an old, red bulge cannot be very large. 



In conclusion, since bulges in Sbc/Sc spirals seem to 
display, at odds with earlier types, a continuity of properties 
with the inner regions of discs (Carollo et al. 1999); the true 
global M/L ratio of an Sbc/Sc spiral should lie somewhere 
between the alternative values given in Table 6 and in Fig 25. 



8.1 M/L offsets with Hubble type: the bulge 
contribution 

In this section we recompute the offsets in stellar M/L ratios 
- and correspondingly in the TF relation - as a function of 
Hubble type, discussed in §5.1, including the effect of bulges. 

The ^-parameter we used as Hubble type indicator, is 
an indicator of the SFH in the disc (see §2) , hence the offsets 
in stellar M/L ratios discussed in §5.1 are essentially due to 
the differences in disc properties. However, the bulge also 
contributes to the global M/L ratio, possibly in a different 
fashion in different Hubble types. It is not clear beforehand 
how the bulge contribution might change the disc-based off- 
sets in §5.1. On one hand, the B/V luminosity ratio increases 
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(1) 

Hubble type 


(2) 
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Sa/Sab 

Sb 
Sbc/Sc 


0.1 
0.35 
0.9 


Salpeter 
Salpeter 
Salpeter 


2.70 
2.09 
1.52 


3.77 
3.77 
3.77 


0.45 
0.35 
0.25 


0.63 

0.63 
0.62 


3.03 
2.53 
1.97 


1.54 
1.28 

1 


0.47 
0.27 



0.99 
0.80 
0.61 


1.31 
1.31 
1.31 


0.48 
0.38 
0.29 


Sa/Sab 

Sb 
Sbc/Sc 


0.1 
0.35 
0.9 


Kroupa 
Kroupa 
Kroupa 


1.76 
1.44 
1.11 


2.35 
2.35 
2.35 


0.45 
0.35 
0.25 


0.60 
0.57 
0.53 


1.94 
1.68 
1.36 


1.43 
1.23 
1 


0.39 
0.23 



0.64 
0.55 
0.45 


0.82 
0.82 
0.82 


0.47 
0.38 
0.29 


Sa/Sab 

Sb 
Sbc/Sc 


0.1 
0.35 
0.9 


Kcnnicutt 
Kcnnicutt 
Kcnnicutt 


1.29 
1.02 
0.76 


1.80 
1.80 
1.80 


0.45 
0.35 
0.25 


0.63 
0.62 
0.59 


1.45 
1.22 
0.97 


1.41 
1.26 
1 


0.37 
0.25 



0.47 
0.39 
0.31 


0.61 
0.61 
0.61 


0.48 
0.39 
0.30 


Sa/Sab 

Sb 
Sbc/Sc 


0.1 
0.35 
0.9 


Larson(l) 
Larson(l) 
Larson(l) 


1.52 
1.13 
0.79 


2.16 
2.16 
2.16 


0.45 
0.35 
0.25 


0.64 
0.67 
0.69 


1.72 
1.40 
1.06 


1.62 
1.31 
1 


0.52 
0.29 



0.55 
0.43 
0.32 


0.74 
0.74 
0.74 


0.48 
0.39 
0.29 


Sa/Sab 

Sb 
Sbc/Sc 


0.1 
0.35 
0.9 


Larson(2) 
Larson(2) 
Larson(2) 


1.56 
1.22 
0.91 


2.18 
2.18 
2.18 


0.45 
0.35 
0.25 


0.63 
0.62 
0.60 


1.75 
1.47 
1.17 


1.50 
1.26 
1 


0.44 
0.25 



0.57 
0.47 
0.37 


0.75 
0.75 
0.75 


0.47 
0.38 
0.29 


Sa/Sab 

Sb 
Sbc/Sc 


0.1 
0.35 
0.9 


Chabrier 
Chabrier 
Chabrier 


1.72 
1.29 
0.91 


2.54 
2.54 
2.54 


0.45 
0.35 
0.25 


0.67 
0.69 
0.70 


1.97 
1.61 
1.24 


1.59 
1.30 
1 


0.51 
0.29 



0.62 
0.49 
0.36 


0.87 
0.87 
0.87 


0.48 
0.39 
0.29 



Table 7. Total M/L ratios, M/L and magnitude offsets as a function of Hubble type, including the contribution of the bulge. (1) Hubble 
type. (2) Typical b- parameter for the disc. (3) Adopted IMF (mass limits as in §3). (4) M/L ratio in the I-band for the disc, from 
the simple models in §5. (5) M/L ratio in the I— band for the bulge, considered as a SSP of solar metallicity and age 10 Gyr. (6) B/T> 
luminosity ratio in the I— band. (7) Actual B/T> mass ratio. (8) Global M/L ratio in the I-band for the galaxy. (9) Offset in M/L ratio 
with respect to Sbc/Sc galaxies. (10) Offset in magnitude with respect to Sbc/Sc galaxies. (11) M/L ratio in the K-band for the disc 
(from the simple models in §5). (12) M/L ratio in the K-band for the bulge (considered as a SSP of solar metallicity and age 10 Gyr). (13) 
B/V luminosity ratio in the K-band corresponding to the mass ratio in column 7, or cquivalcntly to the adopted (B/V)j in column 6. 



toward earlier Hubble types; hence one might expect that 
bulges in earlier Hubble types are more prominent and con- 
tribute more to increasing the global M/L ratio of galaxies; 
this would enhance the systematic offset in M/L of earlier 
types with respect to Sbc/Sc. On the other hand, in earlier 
Hubble types the stellar populations in the disc per se are 
on average older and redder, so that the difference in M/L 
ratio between disc and bulge is less significant; hence one 
might expect that including the bulge has a lesser impact 
on the global M/L ratio for early types than for late types, 
reducing the net offsets. 

In Table 7 we estimate the offsets in global M/L ra- 
tios as a function of Hubble type, adding the effect of the 
bulge to the simple disc models computed in §5. The M/L 
ratios for the discs (column 4) are taken from §5, those of 
the bulges (column 5) are derived from 10 Gyr old SSPs of 
solar metallicity. The B/V ratios in the I-band as a func- 
tion of Hubble type (column 6) are from Sommer-Larsen 
et al. (2003). Computing the actual mass ratios (column 7) 
shows that the B/V ratio in mass does not depend strongly 
on Hubble type, hence the trends in B/V luminosity ratios 
are mostly due to different stellar populations and SFHs in 
the discs — as suggested by Kennicutt et al. (1994). For 
all Hubble types (and all IMFs) the typical B/V mass ra- 
tio is around 60%. However, we remind that this estimates 
are derived considering the bulges of all Hubble types to be 
similar: 10 Gyr old red objects. As discussed above, this is 
a good estimate for Sa to Sb galaxies, less so for Sbc/Sc 
galaxies; if the latter have younger bulges, the mass B/V 
ratio may still have a real trend with Hubble type. Besides, 



what we are discussing here are gross average trends, but 
the scatter in B/V ratios is known to be large (de Jong 
1996a; Sommer-Larsen et al. 2003 and references therein). 
In particular, bulge-less disc galaxies are also observed, in 
particular among late Hubble types. 



Including the effect of the bulge, the typical stellar M/L 
ratio is ~1.5 times larger in Sa/Sab galaxies than in Sbc/Sc 
galaxies (column 9). The corresponding systematic offsets in 
magnitude with respect to the Sbc/Sc zero point are about 
0.25-0.3 mag for Sb galaxies and 0.4-0.5 mag for Sa/Sab 
(column 10). In summary, the effect of the bulge is to reduce 
the offsets expected from the sole differences in disc SFH 
(§5). The offsets however remain significant. 



The present estimate is a lower limit to the expected 
offsets, since the effect of bulges is maximized. As dis- 
cussed in §8, the I-band B/V ratios from Sommer-Larsen 
et al. (2003) are likely to be overestimated, because derived 
from a decomposition performed with ar' profile — hence 
in Table 7 we are overestimating the effects of bulges. In ad- 
dition, if bulges in Sbc/Sc galaxies are younger, the global 
M/L ratios of these galaxies are lower than listed in Table 7, 
and the difference to earlier types becomes more significant. 
All in all, the typical magnitude offsets should lie between 
the values in Table 7 and the estimates given in §5 from the 
pure disc contributions. 
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8.2 The bulge— to— disc ratio in the K— band 

It is very common nowadays to assess issues like the stellar 
mass distribution in galaxies and the B/T> ratio from obser- 
vations in NIR bands, supposedly better mass tracers than 
blue or optical bands, since an important contribution to the 
light in these bands comes from the underlying old stellar 
component. We stressed in §5.2 the uncertainties in estimat- 
ing absolute M/L ratios in NIR bands, especially due to the 
role of AGB stars, a complex phase with uncertain modelling 
in population synthesis. However, as an estimate of relative 
mass content, NIR bands are still useful. For instance, in 
Fig. 8 different treatments of the AGB phase (G-SSPs vs. 
M-SSPs) seem to result in overall offsets of the K-band M/L 
ratios, while the trend as a function of the 6-parameter is 
preserved. Therefore comparing the NIR luminosity of two 
objects, we still get an indication of their stellar mass con- 
tent in relative terms, while deriving their absolute stellar 
mass from the M/L ratio would be more model-dependent. 

Since a lot of focus is set nowadays on NIR observations 
and surveys (e.g. 2MASS, DENIS), it is useful to present 
here the B/T> ratios we predict in the K-band, and discuss 
to what extent the K-band luminosity can be considered a 
direct tracer of stellar mass. 

In Table 6, we list the K-band M/L ratios for our model 
Sbc/Sc discs (column 8) and for the bulges, modelled as 
10 Gyr old SSPs of solar metallicity (column 9). We verified 
that changing the typical metallicity of the bulge by a factor 
of two with respect to solar, would change the bulge M/Lk 
ratio by less than 10%. In column 10 we list the resulting 
B/T> ratios in the K-band when the mass ratios from col- 
umn 6 arc assumed; these are the (B/T>) K ratios correspond- 
ing to the adopted (B/T>) 1 =25% in the I-band, for Sbc/Sc 
galaxies. In all cases, the typical (B/T>) K =30%, which is in 
good agreement, within the large scatter, with the findings 
of de Jong (1996a) for Sbc/Sc galaxies in his r 3 bulge-disc 
decompositions in the K-band. Still, the K-band B/T> ra- 
tio is not a direct indicator of the corresponding mass ratio, 
due to the typical differences in stellar populations and M/L 
ratio between bulge and disc (cf. column 10 vs. column 6). 

In Table 7, columns 11 to 13, we list the K-band M/L 
ratios for discs of different Hubble type (from the simple 
models in §5), the M/L ratios of bulges, and the K-band 
B/T> ratios corresponding to the mass ratios in column 7, 
i.e. corresponding to the [B/T>) 1 ratios adopted in column 6. 
The typical B/T> ratio in the K-band is 30% for late-type 
spirals, 40% for Sb spirals, and 50% for early type, Sa/Sab 
spirals. These values and trends with Hubble type are in ex- 
cellent agreement with de Jong (1996a) for his r* bulge-disc 
decompositions; however, we are speaking of average trends 
but the observational scatter is huge (de Jong 1996a). For 
no Hubble type, the K-band B/T> ratio can be considered 
a direct indicator of the corresponding mass ratio (cf. col- 
umn 13 vs. column 7), so that to derive the B/T> mass ratio 
one should always account for the different stellar popula- 
tions and M/L ratios between bulge and disc. 

9 SUMMARY AND CONCLUSIONS 

We reviewed in the introduction some theoretical and ob- 
servational evidence that the M/L ratio of the stellar com- 



ponent in late-type spirals should be low: M/Li <1 in the 
I-band. This paper is devoted to discussing whether M/L 
ratios so low are compatible with current understanding 
of the stellar populations in galaxies. We explained in §2 
why we should concentrate on late type spirals like Sbc/Sc, 
characterized by a SFH with "birthrate parameter" b >0.8, 
and we argued the the I-band is the optimal band for esti- 
mates of M/L ratios. In this study we considered the "clas- 
sic" Salpeter IMF, mostly for the sake of comparison, and 
other more "bottom-light" IMFs favoured by independent 
studies: the Kroupa IMF, the Kennicutt IMF, the Larson 
IMF (also in a modified version) and the Chabrier IMF (§3). 
We calculated up-to-date sets of SSPs for a wide range of 
metallicities with these different IMFs, to compute the cor- 
responding galactic photometry (§4). A simple approach, 
based on one-zone models with exponential SFHs and fixed 
metallicity, was presented in §5 showing that, aside from the 
Salpeter IMF which has a high M/L ratio, all the other IMFs 
imply M/Li <1 for Sbc/Sc galaxies. Our models also pre- 
dict a significant difference in typical M/L ratios — hence 
in the zero-point of the TF-relation — between spirals of 
different Hubble type. We used our simple models to show 
that the predicted M/L ratio in the I-band is quite robust 
with respect to uncertainties in the current SSPs (mostly the 
contribution of the AGB phase), in contrast to NIR bands. 
Furthermore, we assessed uncertainties related to spectral li- 
braries. We estimated our M/L ratios to be robust within a 
10% or better for the typical metallicities of spiral galaxies, 
and to be conservatively high. 

For a more realistic analysis, we developed chemo- 
photometric models with gas infall and radial gradients in 
SFH and metallicity, to discuss M/L ratios in connection 
with the chemical properties of spiral discs of Sbc/Sc type. 
In the overall the M/L ratios from full chemo-photometric 
models agree quite well with the predictions of simple expo- 
nential models, but some further insight is obtained when 
considering the corresponding chemical evolution. 

Due to the low locked-up fractions, some of the bottom- 
light IMFs (Kennicutt, Larson and Chabrier) result in high 
metal yields, which are hardly compatible with the observed 
metallicities and gas fractions — unless metallicities are 
heavily underestimated due to dust depletion, but this seems 
unlikely (§7.5). Either we reduce the typical metal produc- 
tion (by means of suitable tuning of the upper mass limits 
of the IMF, or by allowing for fallback of metals onto a 
black hole after the SN explosion), or we need to invoke 
dispersal of metals from disc galaxies into the intergalactic 
medium, a behaviour reminiscent of cluster galaxies enrich- 
ing the intra-cluster medium. In any case, the resulting M/L 
ratio for Sbc/Sc discs is M/L/ ~0.7-0.9. 

In other cases (Kroupa and modified-Larson IMFs), the 
slope of the IMF at high masses is quite steep and the metal 
production turns out to bee low; then the net yield must be 
increased to obtain Sbc/Sc discs with the observed metal- 
licities. This is accomplished by tuning the lower mass end 
to larger values, thereby decreasing the locked-up fraction 
but also, correspondingly, the M/L ratio. Also in this case, 
we find M/Li ~0.8 for Sbc-Sc discs. 

Although uncertainties in the gas fractions (a basic con- 
straint for chemical evolution models) and in the theoretical 
stellar yields may hamper the quantitative details of the 
models, e.g. the exact tuning on the IMF mass limits, the 
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predicted M /L ratios appear to be rather robust with respect 
to these details. We consistently find that low M/L ratios, 
M/Li ^0.8 for the typical SFHs of Sbc/Sc discs, naturally 
follow from currently popular, "bottom-light" models of the 
stellar IMF. Even after including the contribution of bulges 
to the overall M/L ratio of disc galaxies, we globally find a 
typical M/L/ < 1 with the "bottom-light" IMFs. 
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APPENDIX A: THE ZERO POINT OF THE 
I-BAND TF RELATION 

One of the motivations for this study was the need for a 
low M/L ratio in galactic discs to match the observed TF 
relation (§1). However, there are differences in the zero point 
of the TF relation among different samples. 

In Fig. 26 we compare different I-band TF relations 
based on large samples of cluster and/or field galaxies, as 
rescaled by their authors to the respective kinematic zero- 
points. The TF relation by Giovanelli et al. is based on a 
sample of 24 clusters, and the zero-point is set by requiring 
that the galaxies in the 14 most distant ones move on aver- 
age with the Hubble flow (negligible peculiar velocities). A 
better kinematic zero-point, based on a sample of 52 distant 
clusters, has been determined by Dale et al. (1999), making 
the TF relation 0.1 mag dimmer than the one by Giovanelli 
et al. Still, even with this correction there is a mismatch 
with other samples — see Fig. 26, top panel. 

The sample by Han & Mould (1992 and companion pa- 
pers) as reanalyzed in the MARK III catalogue by Willick 
et al. (1995), is also a large deep sample of 24 clusters, with 
the zero point set by requiring that the global sample is 
on average moving as a pure Hubble flow. The sample by 
Mathewson, Ford & Buchhorn (1992), also reanalyzed in the 
MARK III catalogue, is a large sample of relatively nearby, 
mostly field spirals, normalized to the Han & Mould sample 
in line widths and kinematic zero point (Willick et al. 1996) ; 
compared to the other two TF relations, this is more influ- 
enced by peculiar motions and bulk flows due to the vicinity 
of the galaxies. The TF relation by Giovanelli et al., which 
we used as a reference in Fig. 1, is the brightest. The large 
offset between the sample by Giovanelli et al. and that by 
Han & Mould can be explained partly by the different defi- 
nitions of line width used: at 50% and 20% of the maximum 
HI flux, respectively (W 50 and W20; Willick 1990; Han & 
Mould 1992). The two estimates can differ by up to a 10% 
(Aaronson et al. 1986; Macri et al. 2000; Sakai et al. 2000), 
or 0.04 in logarithm. Since the TF relation is so steep, mi- 
nor systematic offsets in the estimate of the circular velocity 
become crucial for the zero point in luminosity. 

These offsets in velocity widths are unimportant when 
the TF relation is used as a relative distance indicator or 
to map the peculiar velocity field (provided the definition is 
consistent within the sample), but they correspond to large 
effects in luminosity when the TF relation is used to test 
galactic models. The question obviously arises, which TF 
relation should galactic models be compared to. In Sommer- 
Larsen & Dolgov (2001), Sommer-Larsen et al. (2003) the 
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Figure 26. Comparison of different I— band TF relations with 
kinematic zero— points, from Giovanelli et al. (1997b), from the 
MARK III catalogue (Willick et al. 1995, 1996) and from 
Courteau et al. (2000). Middle panel: relations rescaled to a com- 
mon W50 (see text). Bottom panel: inverse MARK III TF rela- 
tions, rescaled to W50 ( see text). 

circular velocity of the simulated disc galaxies is defined as 
the circular velocity at R = 2.2 Rd- For Milky Way sized 
galaxies this corresponds quite well to the maximum circu- 
lar velocity, as obtained from optical rotation curves (Sack- 
ett 1997; Courteau 1997). As to radio line widths, W 50 
is the one that best matches twice the optical maximum 
velocity (Buchhorn 1992; Mathewson et al. 1992; Courteau 
1997), and it is used in the Giovanelli et al. sample — with 
some further corrections to recover the true maximum rota- 
tional velocity (Giovanelli et al. 1997a). In the middle panel 
of Fig. 26 we rescale the two TF relations of the MARK 
III samples from W20 to W50, dividing by 0.91 (Macri 
et al. 2000) or applying the corrections by Willick (1991), 
as in Han & Mould (1992). These samples are now closer in 
zero point to the Giovanelli sample, yet with different con- 
versions the difference in zero-point can reach 0.4 mag for 
Milky Way sized galaxies (log W ~2.6), with the Giovanelli 
et al. TF being the brightest for V c > 130 km/sec. 

Other effects may contribute to these magnitude offsets, 
most notably different internal dust corrections (Giovanelli 
et al. 1995 vs. Willick et al. 1996) or a different mixture of 
Hubble types. The TF relation of Giovanelli et al. has been 
corrected to apply to late-type spirals, while the other re- 
lations are averaged over Hubble type. All samples consist 
mostly of galaxies between Sb and Sc, but the systematic 
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Figure 27. Comparison of different I— band TF relations with 
locally calibrated zero-points; from Giovanelli et al. 1997b; Tully 
et al. 1998; Sakai et al. 2000; Mathewson ct al. (1992); Buchhorn 
(1992). 



offset in magnitude of Sb galaxies may be larger than esti- 
mated by Giovanelli et al. (see our model predictions in §5.1 
and §8.1). 

On the other hand, the Giovanelli et al. TF relation 
agrees well with the one by Courteau et al. (2000) , based on 
~300 galaxies at a distance of about 6000 km sec -1 (shown 
only in the top panel of Fig. 26 for the sake of clarity). 

Fig. 27 compares available I-band TF relations in ab- 
solute magnitude, calibrated using galaxies with Cepheid- 
based, absolute distances. All relations are expressed as a 
function of the maximum rotational velocity Vmax- 

Here the Giovanelli et al. TF relation is plotted accord- 
ing to their calibration using a local sample with known 
Cepheid distances, corresponding to h — 0.69 (Giovanelli 
et al. 1997b), hence very close to the relation plotted in 
Fig. 1 where h = 0.65 was used. Their linewidths W50 are 
corrected so as to give directly a good estimate of twice Vmax 
(Giovanelli et al. 1997a). 

The Tully et al. (1998) relation is based on two clusters 
(Ursa Major and Pisces) and calibrated with local galaxies of 
known Cepheid distances. To convert their linewidth Wr to 
Vmax we adopted the corrections by Tully & Fouque (1985). 

The Sakai et al. (2000) TF relation is based on a com- 
bination of a local sample (with Cepheid distances obtained 
from the HST Ho Key Project) and the sample by Giovanelli 
et al. We assume here that the linewidth they use, W50, well 
corresponds to twice Vmax — see above. 

The Mathewson TF relations are those defined by 
Mathewson et al. (1992) for Fornax, and by Buchhorn 
(1992) combining clusters and groups in the Mathewson 
sample (hence these TF relations differ from the global 
Mathewson/MARK III TF in Fig. 26). We rescaled both 
relations to the recent, Cepheid-based estimates of the For- 
nax distance modulus, 31.5 on average (Mould et al. 2000). 
Their TF relations are expressed directly in terms of Vmax, 
estimated mostly from optical rotation curves. 

In Fig. 27 also, the Giovanelli et al. TF relation is the 
brightest on average. There is very good agreement with the 
Tully et al. relation, but large differences in slope and/or 



zero point (~0.6 mag) exist in comparison to the other sam- 
ples. 

Part of the differences in zero-point are related to the 
adopted technique for data analysis. One can derive "di- 
rect" TF relations (least square fits minimizing the scatter 
in magnitude, with logarithmic velocity as independent vari- 
able), "inverse" TF relations (least square fits with magni- 
tude as independent variable), or bivariate TF relations (as 
in Giovanelli et al. 1997a). For the same sample, the direct 
TF relation tends to have a shallower slope, hence be dim- 
mer at the bright end, than the inverse TF relation (e.g. 
Willick et al. 1995, 1996). Fig. 26 and Fig. 27 seem to reveal 
good agreement between bivariate (Giovanelli) and inverse 
(Courteau, Tully) TF relations drawn from different sam- 
ples. The Sakai TF relation is also bivariate; it is steeper 
than the Giovanelli et al. but the luminosities are compara- 
ble at the high velocity end. Direct TF relations tend instead 
to be dimmer (MARK III in Fig. 26, top and middle panels; 
Mathewson in Fig. 27). For the MARK III samples, inverse 
TF relations are also available (Willick et al. 1995, 1996). In 
the bottom panel of Fig. 27 we plot the inverse MARK III 
TF relations, also rescaled to a common W50 with the Gio- 
vanelli TF relation. There is now better agreement between 
the luminosities at the bright end. 

Various arguments favour an inverse or bivariate fit over 
a direct fit (Schechter 1980; Tully 1988; Giovanelli 1997a). 
In this sense we may disregard the dimmer direct TF re- 
lations and consider the Giovanelli et al. as in reasonable 
agreement, within uncertainties, with other (bivariate or in- 
verse) TF relations. 



APPENDIX B : COLOUR CORRECTIONS TO 
THE STELLAR M/L RATIO 

In this appendix we discuss the relation between the M/L 
ratio of the stellar (+remnants) component of a galaxy, and 
its colours — hereinafter MLC relation. MLC relations are 
useful, for instance, to colour-correct the TF relation to 
a common zero point for all Hubble types, basing on the 
colours of the individual galaxies (Kannappan et al. 2002), 
or to estimate the radial variation of the stellar M/L ratio 
in a disc, when deriving the mass profile deconvolving the 
baryonic and the dark component (Kranz et al. 2003). The 
luminosities and colours computed here do not include dust 
effects, hence observational colours and magnitudes should 
be dust-corrected before our MLC relations can be applied. 
However, the reddening+dimming effects of dust result in 
a vector that runs almost parallel to the MLC relation, so 
that to a first approximation dust has a minor influence on 
stellar mass estimates from the MLC relation (Bell & de 
Jong 2001, hereinafter BdJ). 

For simplicity, in this appendix we consider only IMFs 
with the "standard" mass limits given in §3: [0.1-100] Mq 
for the Salpeter, Kroupa and Kennicutt IMFs and [0.01- 
100] Mg for the Larson, modified Larson and Chabrier IMFs 
(in these cases whether the lower limit is 0.01 or 0.1 Mq has 
negligible impact on the M/L ratios). Other mass ranges, 
considered for example in Table 2 when tuning the level of 
metal enrichment, would correspond to simple changes in the 
zero-point of the MLC relation — as will become apparent 
from what follows. 
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Fig. 28 shows some examples of MLC relations for three 
representative IMFs among those considered in this paper. 
The plot shows results from our chemo-photometric models 
of disc galaxies, both the calibrated models in Table 2 and 
several other models computed during the calibration pro- 
cess. Open symbols represent results for the individual one- 
zone annuli that constitute our disc models (see §6); solid 
symbols represent results for multi-zone regions, namely 
they are the integrated luminosities and colours of the sum 
of several annuli (e.g. the "bulge region" within R < 0.5 hs, 
or the "disc region" R > 0.5 /is, or the global galaxy, etc.; 
see §6.3). Asterisks represent results for the "optical disc 
region" (0.5 Kb < R < 3 hs) discussed extensively in the 
paper (§6.3). Henceforth, points falling in the same colour- 
M/L ratio region in the plots result from a variety of star 
formation and chemical enrichment histories, both one-zone 
and composite. This is important to see if colours can be 
good tracers of the underlying M/L ratio variations, with 
no need of independent knowledge of the star formation and 
metal enrichment history. 

A well defined linear relation between colour and log- 
arithm of the M/L ratio is apparent in the plots (see also 
BdJ), at least down to a certain limit in colour (e.g. ~0.5 
in B — V); blueward of this, models deviate toward M/L ra- 
tios lower (higher in the K band) than the linear fit. In the 
"linear range", the MLC relation is rather tight in spite of 
the variety of SFH and chemical evolution history that can 
result in a given colour. The worst correlations are obtained 
with V — K, and in general with colours involving the K 
band (B — K, R — K etc.) or in MLC relations involving the 
K luminosity. The reason is that in optical colours, the age- 
metallicity degeneracy favours convergence into a tight MLC 
relation, while in IR bands and colours, age and metallicity 
effects are less degenerate (see BdJ, their Fig. 2 and related 
comments). Hence the K band (and NIR bands in general) 
offers the advantage of a smaller variation of the M/L ratio 
with colour, Hubble type, SFH and metal enrichment his- 
tory, but colour corrections and MLC relations based on the 
K band are less tight and compelling. 

Fig. 28 shows results for the Salpeter, Kroupa and Ken- 
nicutt IMFs as representative cases. To a first approxima- 
tion, the slope of the MLC relation does not depend on the 
IMF, which just sets the zero-point; here we confirm the 
findings of BdJ. For all the 6 IMFs considered in this paper 
we did verify that, when the linear fit is performed over the 
same colour range, the slope of the MLC relation is indepen- 
dent of the IMF within the uncertainty given by the scatter. 
We provide in Table 8 the typical, IMF-independent slopes 
of the MLC relations, as well as the zero-points correspond- 
ing to the different IMFs. We express the zero-points as the 
logarithmic M/L ratios corresponding to certain values of 
the colour; specifically: 

(B — V) : log (J^j = Sx [( B -V)- 0.6] + qx 
(B-R) : log(^;) = sx[(B - R) - 1.0] + q x 
(V-I) : log (il) = s x [(V -I)- 0.9] + qx 
(V-K) : log(^=sx[(V-K)-2.5}+q x 





colour 


range 


SB 


sv s R 


SI 


SK 




B-V 


0.55-0.90 


1.69 


1.29 1.11 


0.97 


0.73 




B-R 


0.95-1.45 


1.17 


0.89 0.77 


0.68 


0.51 




V-I 


0.88-1.14 


2.17 


1.66 1.43 


1.26 


0.96 




V-K 


2.4-3.0 


1.23 


0.97 0.84 


0.76 


0.57 




colour 


QB 


IV 




11 








Salpeter 


IMF [0.1-100] M Q 






B 


- V 


0.311 


0.341 


0.317 


0.260 


-0.115 


B 


- R 


0.272 


0.312 


0.292 


0.236 


-0.134 


V 


- I 


0.246 


0.292 


0.275 


0.220 


-0.147 


V 


- K 


0.093 


0.169 


0.168 


0.123 


-0.222 






Kroupa IMF [0.1-100] M 






B 


- V 


0.089 


0.120 


0.099 


0.045 


-0.298 


B 


- R 


0.061 


0.099 


0.080 


0.027 


-0.311 


V 


- I 


0.068 


0.104 


0.086 


0.032 


-0.309 


V 


- K 


0.002 


0.048 


0.036 


-0.015 


-0.343 



Kennicutt IMF [0.1-100] M 



B - 


V 


-0.002 


0.028 


0.007 


-0.050 


-0.421 


B - 


R 


-0.034 


0.004 


-0.014 


-0.068 


-0.435 


V - 


I 


-0.041 


-0.002 


-0.020 


-0.074 


-0.438 


V - 


K < 2. 


S -0.182 


-0.114 


-0.117 


-0.162 


-0.505 






Larson IMF [0.01-100] M Q 






B - 


V 


0.081 


0.111 


0.088 


0.029 


-0.359 


B - 


R 


0.045 


0.085 


0.065 


0.010 


-0.373 


V - 


I 


0.023 


0.068 


0.050 


-0.004 


-0.383 


V - 


K < 2.. 


S -0.155 


-0.075 


-0.075 


-0.118 


-0.466 






modified Larson IMF [0.01-100] 


M Q 




B - 


V 


0.064 


0.095 


0.073 


0.016 


-0.356 


B - 


R 


0.031 


0.069 


0.051 


-0.004 


-0.371 


V - 


I 


0.025 


0.065 


0.047 


-0.006 


-0.374 


V - 


K 


-0.123 


-0.055 


-0.059 


-0.104 


-0.446 






Chabrier IMF [0.01 


-100] M Q 




B - 


V 


0.131 


0.161 


0.138 


0.080 


-0.303 


B - 


R 


0.095 


0.134 


0.115 


0.060 


-0.318 


V - 


I 


0.078 


0.121 


0.103 


0.049 


-0.325 


V - 


K < 2. 


S -0.107 


-0.028 


-0.029 


-0.072 


-0.419 



Table 8. Coefficients for the MLC relations. The slopes sx are 
independent of the IMF, while the zero-points qx depend on the 
IMF. "Range" indicates the colour range where the linear fit is a 
good representation of the MLC relation; the lower limit is set by 
the "deviations" from a linear relation (see text and Fig. 28), the 
upper limit by the reddest colour probed by the models. For the 
Kennicutt, Larson and Chabrier IMFs, more strict upper limits 
apply in V-K, as indicated (see text). 



Table 8 lists the slope s and zero-point q for the MLC re- 
lations in the B,V,R,I,K bands; also given are the color 
ranges over which the linear fits are valid. MLC relations 
with a different colour base (e.g. V — R,B — I , etc.) can be 
obtained as a linear combination of these. 

The models with the Kennicutt IMF in Fig. 28 do not 
reach so red colours as the other two sets, because a signifi- 
cant SFR at late times is maintained by the large returned 
gas fraction from the stellar populations (1 — a — 40 — 50% 
for the Kennicutt, Larson and Chabrier IMFs, Table 3). 
As a consequence, low 6-parameters (corresponding to red 
colours) are hardly achieved with these IMFs, unless gas 
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Figure 28. MLC relations from chemo-photometric models of disc galaxies with three representative IMFs. Open symbols: one-zone 
models for individual disc annuli; solid symbols: multi-zone models (integrate photometry of bulge region, disc region, or global galaxy); 
asterisks: "optical disc" region (see text). The solid, short dashed and long dashed lines represent the linear fits for the different IMFs, 
over the respective relevant colour range. The thin dotted line is the linear fit from Bell & de Jong (2001). 



outflows occur, not included in our disc models. The slopes 
sx given in Table 8 have been derived from models with 
the Salpeter, Kroupa and modified-Larson IMFs, for these 
models extend to redder colours so the slopes could rely on 
a wider colour range. 

For the Kennicutt (and Larson and Chabrier) case, the 
linear relation seems to break down also toward the red in 
V — K (Fig. 28, lower right panel), as well as in other colours 
involving the K-band (B — K, R — K etc.). It is not clear 
if this is due to a selection effect — all models with these 
IMFs are blue in optical colours — or whether it reflects the 
fact that the evolution of V — K is more sensitive to the IMF 
than bluer colours (Vazdekis et al. 1996, 1997). Anyways, in 
such cases we give both lower and upper limits to the "linear 
range" where the linear fit can be applied. 

A thorough discussion of MLC relations and their im- 
plications on the TF relation is given by BdJ. The dotted 
lines in Fig. 28 are the linear MLC relations derived by BdJ 
for their favourite set of galactic models (their Table 1). 
Overall, their slopes are in good agreement with ours, albeit 
slightly steeper (with the exception of the MLC relations 



in V — K). Our slopes are somewhat intermediate between 
those of their full galactic models (their Table 1 and Fig. 28) 
and those of their simple exponential models at fixed metal- 
licity (their Table 4). We impute this slight difference to 
a tighter calibration in metallicity for our galactic models 
(§6.3), while BdJ select their models mostly on the base of 
their K-band magnitude and central surface brightness. 

Their normalization of the MLC relation, based on max- 
imal disc arguments, corresponds roughly to a Kroupa IMF. 
Other IMFs with lower M/L ratios, such as the Kennicutt 
IMF (Fig. 28), would then imply slightly submaximal discs. 
The Salpeter IMF seems to be excluded from maximal disc 
arguments (BdJ). 

The deviation from the MLC relation at blue colours 
seen in Fig. 28 is not apparent in the results by BdJ. How- 
ever, in their Fig. 10 their galactic models with a more reg- 
ular SFH (their closed box, infall, and outflow case), which 
are probably more similar to our chemo-photometric mod- 
els, do show a slight bend-over in the MLC relation and do 
not stretch to very blue colours. In their favourite set of 
models, small bursts of star formation are included and the 



© 2003 RAS, MNRAS 000, 000-000 



32 Portinari et al. 



bend-over is wiped out by the larger scatter. The different 
behaviour toward blue colours is probably due to the fact 
that, while their models are global galactic models selected 
so as to have K band magnitudes and central surface bright- 
nesses representative of spiral galaxies, our blue models (say, 
bluer than B-R~0.9 for instance) corresponds to individual 
external annuli of our galactic disc models, characterized by 
overall young ages and low metallicities. 

Our steeper slope in the MLC relations in V-K is prob- 
ably due to the fact that it is defined over a narrower 
colour range, out of which our MLC relation tends to flatten 
(Fig. 28, bottom right panel). 

Considering that both our galactic chemical evolution 
models and our population synthesis code differ from BdJ, 
we regard our results as in good agreement with theirs, 
which supports the use of MLC relations as a robust theo- 
retical prediction. 
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